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ABSTRACT 


A theory  of  estimation  of  angular  position  and  other  at- 
tributes of  optically  rough  and  smooth  targets  with  a monopulse 
laser  radar  is  developed.  It  applies  to  systems  deriving  in- 
formation about  the  target  position  by  sensing  its  image  with 
an  array  of  noncoherent  detectors.  The  theory  develops  quan- 
titative formulations  of  the  fundamental  limitations  imposed 
upon  measurement  accuracy  by  the  shot  noise  arising  from  both 
the  target  return  and  the  background  radiation,  by  the  detec- 
tor dark  current,  and  by  the  random  fluctuations  of  the  target 


cross  section. 
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I.  INTRODUCTION 

The  potential  of  considerably  higher  accuracy  than  that 
of  microwave  radar  motivates  interest  in  employing  laser  radar 
for  angle  tracking.  With  the  recent  advent  of  high-power  la- 
sers, angle -tracking  laser  radars  are  particularly  attractive 
for  space-borne  applications  where,  due  to  the  absence  of  at- 
mospheric turbulence,  the  potential  for  high  accuracy,  which 
results  primarily  from  employing  short  wavelengths,  can  be  used 
to  full  advantage.  High-accuracy  tracking  used  for  precision 
pointing  of  a communication  laser  beam  from  one  satellite  to 
another  distant  satellite  promotes  efficient  use  of  laser  en- 
ergy in  the  communication  link.  It  has  also  been  pointed  out 
that  precise  knowledge  of  the  relative  positions  of  two  space- 
craft performing  a rendezvous  maneuver  is  difficult  and  often 
impractical  to  obtain  from  ground-based  radars;  accurate  laser 
radars  located  on  board  a spacecraft  could  serve  as  accurate 
guidance  aids. 

In  an  angle -tracking  system,  a sequence  of  direct  measure- 
ments of  target  position  is  fed  into  a tracking  filter,  which 
can  produce  refined  estimates  of  past,  present  and  future  tar- 
get position.  Generally,  the  accuracy  of  the  refined  estimate 
is  proportional  to  the  accuracy  of  direct  measurements,  so  that 
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schemes  for  improving  measurement  accuracy  are  potentially  of 
great  significance. 

Relative  to  the  volume  of  research  on  angular  tracking  of 
unresolvable  targets  with  microwave  radars , initiated  during 
World  War  II,  or  on  passive  tracking  of  resolved  targets  re- 
ported more  recently  in  the  infrared  and  visible  regions,  very 
few  papers  have  been  published  bearing  on  the  problem  dealt 
with  herein. 

When  applied  to  tracking  in  the  optical  region,  such  well- 
known  and  established  microwave  radar  tracking  techniques  as 
amplitude  and  phase  monopulse  need  to  be  reassessed  and  ex- 
tended to  cover  at  least  two  new  aspects.  First,  the  classi- 
cal wave  description  used  successfully  in  microwave  radars  is 
not  adequate:  since  the  energy  of  a single  optical  photon  is 

so  much  greater  than  that  of  a microwave  photon,  there  are  far 
fewer  photons  arriving  at  the  receiver  in  any  case  of  practi- 
cal interest.  One  cannot  assume,  even  in  principle,  that  there 
exists  a noise-free  limit  in  which  it  is  possible  to  find  the 
position  of  an  object  exactly.  The  position  estimator  of  an 
optical  object  is  always  subject  to  error  which,  in  the  most 
optimistic  case,  is  quantum-noise  limited.  Second,  the  theory 
of  angular  position  estimators  of  nonresolvable,  point-type  tar 
gets  must  be  extended:  the  fine  resolution  of  optical  imaging 

systems  allows  one  to  distinguish  various  features  on  the  ob- 
ject, and  tracking  of  resolvable  targets  must  be  considered. 
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Furthermore,  microwave-radar  design  literature  offers  little 
or  no  information  for  tracking  targets  illuminated  with  par- 
tially coherent  radiation,  as  might  be  the  case  where  efficient 
multi-line  HF  lasers  are  used,  or  for  the  systems  using  direct 
detection. 

This  thesis  addresses  the  issue  of  target  angular  posi- 
tion (TAP)  estimation  with  a pulsed  laser  radar.  The  TAP  es- 
timates are  obtained  from  the  targets ' images . The  imaging 
system  receives  ligth  from  both  the  target  that  has  been  illu- 
minated with  the  laser  beam  provided  for  tracking  and  from  uni- 
form, noncoherent  background  radiation.  This  light  is  imaged 
onto  the  detector  array,  which  converts  the  photon  flow  into 
an  image  photocounting  distribution  (IPD),  the  electrical  sig- 
nal used  by  the  image  processor.  The  purpose  of  the  image 
processor  is  to  estimate  accurately  the  image  position  for  each 
received  pulse. 

The  thesis  analyzes  several  image  processors  based  on  the 
derived  estimators  of  the  TAP.  Their  accuracy  is  discussed  and 
compared.  In  this  comparison,  the  beam-splitting  ratio,  which 
is  understood  here  as  the  factor  by  which  the  processor  accu- 
racy can  be  improved  above  the  diffraction-limited  resolution 
of  the  imaging  system,  is  used  as  the  figure  of  merit.  In  the 
next  section  the  scope  of  the  thesis  is  outlined.  The  results 
are  summarized  in  the  last  section  of  this  chapter. 


>\-vVvVr  Vy3 
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1.1  Outline 

The  analytical  results  presented  in  this  thesis  have  been 
developed  in  four  major  steps: 

a.  The  semiclassical  photon -counting  theory  was  used  to 
obtain  a model  for  the  Image  Photoelectron-Counting  Profiles 
produced  at  the  array  output  in  response  to  the  irradiance  dis- 
tribution at  the  image  plane.  This  irradiance  was  assumed  to 
be  a sum  of  two  terms,  representing  the  target  image  and  a uni- 
form, noncoherent  background.  A mathematical  model  for  the 
image  irradiance  that  takes  into  account  the  coherence  of  il- 
lumination, target  surface  roughness,  and  the  illuminator-tar- 
get-imaging system  geometry  was  developed.  Mandel's  transfor- 
mation was  used  to  relate  this  model  to  the  statistical  proper- 
ties of  the  image  photoelectron-counting  profile. 

b.  The  faithfulness  with  which  the  image  photocounting 
distribution  provides  knowledge  about  the  image  irradiance  was 
addressed  by  calculating  the  signal-to-noise  ratio  at  the  out- 
put of  the  individual  detectors  in  the  array.  A general  ex- 
pression for  the  SNR  was  obtained  as  a function  of  the  relevant 
target  and  system  parameters. 

c.  For  the  targets  much  smaller  than  the  diffraction- 
limited  resolution  of  the  imaging  system,  their  angular  posi- 
tion was  estimated,  based  on  image  sensing  by  a four-quadrant 
detector.  The  lower  Cramer-Rao  bounds  were  established  for 
the  estimates.  The  maximum  likelihood  estimates  of  the  angular 
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position  were  derived  for  optically  rough  and  smooth  targets 
and  their  performance  compared  with  the  theoretical  limits. 

d.  The  problem  of  estimation  of  angular  position  of  ex- 
tended targets,  well  resolved  by  the  pupil  aperture,  was  ad- 
dressed by  considering  the  properties  of  the  image  centroid 
estimator.  The  mean  and  variance  of  this  estimator  were  ex- 
pressed in  terms  of  the  covariance  matrix  of  the  photoelectrons 
at  the  outputs  of  the  detectors  that  form  the  elements  of  the 
large  array.  The  properties  of  this  matrix  were  discussed  for 
optically  rough  and  smooth  targets.  The  variance  of  the  cen- 
troid estimator  was  obtained  as  a function  of  relevant  target 
and  system  parameters. 

1 . 2 Summary  of  Results 

The  problems  of  TAP  estimation  with  monopulse  laser  ra- 
dars has  been  analyzed  for  systems  that  derive  the  desired  in- 
formation out  of  a single  pulse  from  the  target's  image,  which 
is  sensed  with  an  array  of  noncoherent  detectors. 

In  Chapter  II,  a general  expression  for  the  signal-to- 
noise  ratio  (SNR)  at  the  output  of  an  individual  detector  in 
the  array  is  obtained  as  a function  of  the  coherence  of  the 
illuminating  light,  the  object  surface  roughness,  the  width  of 
the  telescope  point  spread  function,  and  the  aperture  and  in- 
tegration time  of  the  detector.  The  expression  is  evaluated 
for  several  limiting  cases  of  coherence  of  illumination  and  of 
object  surface  roughness.  The  expression  indicates  that,  for 
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increasing  signal  levels,  the  SNR  is  limited  first  by  back- 
ground noise,  then  by  signal  shot  noise,  and  ultimately  by  the 
target  fading.  The  target  fading  manifests  itself  as  the  ran- 
dom fluctuations  of  the  laser  radar  target  cross  section. 

The  fading  contrast  involves  stochastic,  spatial,  and 
temporal  averaging  of  the  second  moment  of  the  image  irradiance 
or  the  degree  of  coherence,  which  propagates  from  the  rough  tar- 
get to  the  image  plane.  The  roughness  enters  by  way  of  the 
joint  characteristic  function  of  the  surface  statistics.  The 
expressions  obtained  for  the  fading  contrast  reveal  quantita- 
tively the  efficiency  of  various  means  for  improving  fading- 
limited  SNR,  by  temporal  averaging,  by  spatial  averaging,  or 
by  increased  bandwidth  of  illumination.  Plots  of  the  SNR 
against  bandwidth  for  various  roughness  scales,  and  against 
signal  strength  for  various  noise  and  fading  levels,  are  pre- 
sented and  it  is  shown  that  the  SNR  can  be  improved  at  the  ex- 
pense of  spatial  resolution. 

In  Chapter  III,  the  angular  position  of  a target  unre- 
solved by  the  imaging  system  is  estimated  by  measuring  the 
location  of  the  target's  photocounting  image  by  means  of  a non- 
coherent four-quadrant  detector.  The  choice  of  estimator  de- 
pends on  the  statistics  of  the  photocounts,  which  reflect  those 
of  the  energy  deposited  on  the  detectors  by  the  laser  field  re- 
turned from  the  target.  For  an  optically  rough  target,  this 
is  a spatio-temporal  complex  Gaussian  process,  producing  photo- 
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electrons  that  obey  a negative  binomial  distribution  with  M 
degrees  of  freedom.  The  limitations  imposed  by  target  fading 
and  by  signal  and  background  shot  and  dark  current  noise  on 
the  accuracy  of  the  estimate  are  investigated.  The  theoret- 
ical limits  given  by  the  Cramer-Rao  lower  bounds  for  the  stan- 
dard deviations  of  the  estimates  are  derived  for  fading  and 
nonfading  targets . For  several  cases  it  is  shown  that  the 
lower  bounds  are  given  by  the  beam-splitting  formula;  i.e., 
that  the  standard  deviations  are  proportional  to  the  diffrac- 
tion-limited angular  resolution  of  the  imaging  telescope  and 
inversely  proportional  to  the  SNR,  where  the  noise  includes 
the  shot  noise  generated  by  the  target  and  background  radia- 
tion and  dark  current,  and  do  not  depend  on  the  target  fading. 

Maximum  likelihood  estimates  are  derived  for  fading  and 
nonfading  targets  and  their  means  and  standard  deviations  are 
evaluated.  All  these  estimates  share  the  following  properties: 

a.  They  depend  on  the  familiar  monopulse  ratio  which, 
for  the  photocounting  problem  is  given  by  the  ratio  of  the 
difference  of  the  total  number  of  photocounts  at  the  detector 
outputs  to  their  sum.  This  normalization  makes  the  angular 
measurement  independent  of  the  target  range  and  significantly 
reduces  the  impact  of  target  fading  on  its  accuracy. 

b.  They  are  biased.  For  nonfluctuating  targets,  the 
bias  decreases  exponentially  with  the  signal  and  background 
noise.  For  fading  targets  the  bias  decreases  monotonically 
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with  the  signal  at  a rate  that  depends  on  the  number  of  de- 
grees of  freedom  in  the  image  radiation. 

c.  The  standard  deviations  of  the  estimates  approach 
the  theoretical  Cramer-Rao  limits  for  unbiased  estimates  when 
the  signals  are  large  enough  to  make  the  MLEs  unbiased. 

The  last  chapter  discusses  TAP  and  target  size  estimation 
of  extended  targets.  The  image  centroid  estimator,  and  that 
for  its  second  moment,  are  obtained  by  processing  the  image 
photocounting  distribution  at  the  output  of  the  multielement 
detector  array  sensor.  The  expressions  for  the  means  and  co- 
variance  matrices  of  the  estimators  are  derived  and  discussed. 

a.  Both  estimators  are  biased;  the  biases  are  inversely 
proportional  to  the  square  of  the  overall  SNR  at  the  sensor 
output. 

b.  The  rms  errors  are  inversely  proportional  to  the  SNR. 

c.  The  biases  and  rms  errors  of  the  estimators  depend 
on  the  detection  matrix,  comprised  of  the  covariances  between 
the  photocounts  at  the  outputs  of  pairs  of  detectors  in  the 
array.  The  detection  matrix  becomes  diagonal,  giving  the  var- 
iances of  the  photocounts  at  the  individual  detector  outputs, 
for  either  non fluctuating  targets  or  noncoherent  illumination, 

or  also  for  fluctuating  targets  generating  Gaussian  image  fields, 
if  the  detectors  sample  them  at  the  Nyquist  rate. 
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d.  For  diagonal  detection  matrices  and  nonfluctuating 
targets,  the  bias  and  rms  error  of  the  centroid  estimator  de- 
crease with  increasing  laser  power;  for  flutuating  targets 
they  cannot  be  decreased  below  the  limit  imposed  by  laser- 
speckle  noise o This  limitation  improves  with  the  number  of 
degrees  of  freedom  in  the  image,  which  depends  on  the  diffrac- 
tion-limited resolution  of  the  imaging  system. 
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II.  SNR  IN  PHOTOCOUNTING  IMAGES  OF  ROUGH  OBJECTS 
IN  PARTIALLY  COHERENT  LIGHT 


Image  irradiance  distributions  from  objects  illuminated 
with  partially  coherent,  quasi-monochromatic  light,  viewed 
against  a spatially  uniform  background  and  received  with  a 
photosensitive  detector  are  analyzed.  A general  expression  for 
the  signal-to-noise  ratio  at  the  detector  output  is  obtained 
as  a function  of  the  coherence  of  the  illuminating  light,  the 
object  surface  roughness,  the  width  of  the  telescope  point 
spread  function,  and  the  aperture  and  integration  time  of  the 
detector.  The  expression  is  evaluated  for  several  limiting 
cases  of  coherence  of  illumination  and  of  object  surface  rough- 
ness . 


Published  in  Appl.  Optics,  15,  2268-2275  (1976)  with  P.  Diament 
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2.1  INTRODUCTION 

Communication  theory  shows  that  the  rate  of  transmission  of 
information  is  limited  by  two  basic  system  properties:  the  bandwidth 
through  which  the  signal  information  passes,  and  the  noise  inherently 
present  in  all  systems.  Analysis  of  the  trade-off  between  signal-to- 
noise  ratio  (SNR)  and  bandwidth  is  an  important  part  of  the  evaluation 
of  communication  systems. 

The  information  conveyed  about  an  object  by  an  imaging  system 
can  be  encoded  both  in  time  and  in  space.  For  systems  employing 
direct  detection,  it  is  the  spatial  correspondence  between  the  irradiances 
at  the  objects  and  at  their  images  that  is  of  primary  interest.  In 
systems  such  as  laser  tracking  radars,  estimates  of  the  position  of 
laser-illuminated  objects  are  derived  from  their  images.  If  the  images 
are  to  be  observed  with  an  array  of  detectors,  the  SNR  analysis  can 


help  decide  whether  the  array  should  consist  of  many  detectors  with 
small  apertures  (high  complexity,  high  resolution)  or  fewer  detectors 
with  larger  apertures  (less  complexity,  lower  resolution). 

The  generality  of  the  system  to  be  considered  is  suggested  in 
Fig.l.  The  object  is  illuminated  with  quasi-monochromatic,  partially 
coherent  light.  The  surface  of  the  object  has  a height  profile  z)^) 
and  a reflection  coefficient  i)^)  that  are  stochastic  processes.  We 
take  z(xQ)  to  be  a stationary,  Gaussian,  zero-mean  process.  The 
widths  of  the  correlation  coefficients  of  the  two  processes,  Pz  and  , 

are  assumed  sufficiently  large  with  respect  to  the  wavelength  to  justify 
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use  of  the  scalar  theory  of  diffraction.  As  a result  of  scattering  from 
the  object  surface,  the  transverse  coherence  originally  present  in  the 
illuminating  field  is  modified,  so  that  the  field  leaving  the  object  becomes 
random  in  time  and  in  space.  The  imaging  telescope  further  modifies 
the  transverse  coherence  of  the  field  at  the  image  plane.  The  total 
image -plane  ir radiance  distribution,  consisting  of  the  sum  of  object  and 
background  irradiance  contributions,  is  sensed  with  a photodetector, 
which  imposes  additional  noise  processes  associated  with  photon  counting 
and  dark  current.  We  seek  the  SNR  at  the  output  for  an  object  illuminated 
with  partially  coherent,  quasi-monochromatic  radiation,  in  the  presence 
of  uniform,  thermal  background  radiation. 

The  imaging  of  partially  coherent  objects  was  first  studied  by 

1 2 3 

Hopkins  and,  more  recently , by  Ichioka  . Enloe  , Lowenthal  and 

4 5 

Arsenault  , and  Dainty  have  analyzed  the  statistical  properties  of  the 

images  of  coherently-illuminated  diffuse  objects.  They  found  that  the 

noise -like,  speckled  structure  of  the  image  has  a mean  square  intensity 

equal  to  the  square  of  the  mean  intensity.  Enloe,  identifying  the  signal 

with  the  expected  value  of  the  fluctuations  of  the  image  irradiance 

distribution  and  the  noise  with  their  standard  deviation,  suggested  that 

the  SNR  can  be  improved  at  the  expense  of  the  spatial  resolution  of  the 

6 7 

imaging  system.  Elbaum,  Greenebaum  and  King  , and  George  and  Jain 
have  studied  the  improvement  in  SNR  when  multi-frequency  illumination 
is  employed  instead  of  nearly-monochromatic  light.  The  effects  of  sur- 
face roughness  on  the  statistical  distribution  of  image  speckle  intensity 


AY.; 
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have  been  measured  by  Fujii  and  Asakura  . 

9 

Goodman  was  the  first  to  discuss  the  properties  of  laser 
speckle  patterns  as  they  affect  the  performance  of  optical  radars. 

Treating  the  field  scattered  from  a coherently  illuminated  object  as  a 
complex  Gaussian  process  over  space,  Goodman  applied  Mandel's  theory 
of  photodetection  to  compute  the  resulting  photocounting  statistics.  This 
paper  combines  Wolf's***  theory  of  propagation  of  coherence  and  Mandel's** 
photodetection  theory  to  obtain  generalized  expressions  for  SNR. 

12 

The  expressions  obtained  extend  the  one  derived  by  Bures  for 
measuring  degeneracy  of  light,  are  not  restricted  to  Gaussian  light  and 
include  background  and  dark  current  noises.  When  specialized  to 
Gaussian  linearly  polarized  light,  they  reveal  explicitly  the  dependence  of 
the  SNR  on  coherence  of  illumination,  object  surface  roughness,  coherent 
point  spread  function,  aperture  of  detector,  and  integration  time.  Although 
interest  is  directed  to  irradiances  at  the  image  plane,  most  of  the 
results  obtained  in  this  paper  can  be  adapted  to  apply  to  irradiances  at  any 
plane. 

In  the  next  section,  the  semiclassical  theory  of  photodetection  is 
employed  to  obtain  a general  formula  for  the  SNR  at  the  output  of  the 
detector.  In  Section  2 .3  the  SNR  is  discussed  for  large  signal  levels, 
where  the  limitations  are  imposed  by  the  random  fading  of  the  field 
across  the  detector.  Several  limiting  cases  of  illumination  coherence 
and  surface  roughness  are  considered.  Section  2 . 4 treats  the  SNR  for 
small  signal  levels,  where  signal-shot-noise  and  background  and  dark 


’</■ 
V /-*■  J . 
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current  noise  must  be  included.  The  tradeoffs  between  SNR  and  spatial 
resolution  are  considered  in  Section  2.5. 

2 . 2 General  Expression  for  Signal-to  Noise  Ratio 

We  seek  to  estimate  significant  features  of  the  irradiance  at  the 
object  from  observations  of  the  number  of  photocounts  at  the  output  of 
a detector  that  senses  the  image  irradiance  distribution.  We  consider  a 
photosensitive  detector  with  an  aperture  A illuminated  by  the  image 
irradiance  I(xj,  t).  In  general  the  total  number  of  counts  n^  observed 
at  the  output  of  the  detector  during  some  integration  time  T includes 
signal  and  background  photoelectrons,  n , n , and  dark  current  electrons 

b D 9 


nD’  ,e  ■ 


nT  ' nS  + nB  + "D  ' 


The  first  two  terms  are  generated  by  the  energy  deposited  on  the  detector 
by  the  electromagnetic  fields  that  arrive  from  the  object  and  from  the 
background.  The  third  term  is  generated  internally  in  the  detector  and 
depends  only  on  the  detector  quality.  We  associate  the  signal,  S,  with 
the  expected  value,  (n8  ) . The  noise  level,  N,  may  be  defined  as  the 
standard  deviation  of  the  counts,  ct^,  , which  measures  the  random  fluctu-? 
ations  at  the  detector  output.  Since  the  three  terms  in  (1)  are  generated 
by  independent  physical  mechanisms  and  may  be  treated  as  mutually 
independent  stochastic  processes,  the  variances  of  the  signal,  background, 
and  dark  current  counts  simply  add,  so  that  the  noise  level  is 

N * °T  ’ (aS  + + a!D)‘A-  <2» 


v*  * 
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The  semiclassical  theory  of  photoemission  developed  by  Mandel^ 
treats  the  number  of  photoelectrons,  n,  resulting  from  the  total 
electromagnetic  field  energy,  W=WT>  deposited  on  the  detector,  as  a 
nonhomogeneous  Poisson  process  with  mean  and  variance  given  in  terms 
of  those  of  W by 

<n>  = a (W)  (3) 

and 

<7n2  = <n>  + a2a2w  . (4) 

The  energy  is  related  to  the  irradiance  by 

W = J*  dt  J d2x  I(x  t)  (5) 

T A 

and  the  coefficient  of  proportionality,  a = ?/hv  , involves  the  quantum 
efficiency  of  the  detector,  §,  and  the  average  frequency  of  the  field,  V; 
h is  Planck's  constant.  The  first  term  in  (4)  is  the  variance  of  the 
photoelectron  count  associated  with  the  photoelectric  effect  and  represents 
the  intrinsic  fluctuations  in  the  detection  process.  The  second  term, 
which  depends  on  the  fluctuations  of  W,  represents  an  excess  fluctuation 
of  counts  due  to  bunching  of  photoelectrons.  The  process  is  thus  doubly 
stochastic. 

When  the  integration  time,  T , is  much  larger  than  the  coherence 
time  of  the  background  radiation,  the  temporal  fluctuations  are  averaged  out 
and  the  bunching  effect  for  the  background  radiation  is  absent;  i.  e.  , 

°B  = ^nB  ^ ' Similarly,  the  dark  count  is  assumed  to  be  distributed 
according  to  Poisson  statistics^,  with  parameter  (n^),  so  that 


'vv 
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Cj)  - (n  d ) • It:  follows  that  the  SNR  can  be  written  as 


SNR  = 


(ns> 


< ns> 


<ns)+a2  + <nN)]‘ 


where  <nN>  = <-B>+<nD>,  and  W is  the  energy  deposited  on  the 
detector  by  the  field  from  only  the  object.  To  analyze  the  relative 
strengths  of  the  three  contributors  to  the  noise,  a modified  version  of 
(6),  using  (3),  is  useful 


-(xb 


<0  . -1/Z 


+ r + — 

<W>2  < 


nc>  ' 


The  second  term  depends  on  the  spatio-temporal  properties  of  I(xj,t) 

modified  by  the  temporal  and  spatial  integration  at  the  detector. 

For  the  special  case  of  an  aperture  and  integration  time  small 

with  respect  to  the  coherence  time,  the  energy  becomes  W = IAT  and 

2 

the  middle  ratio  can  be  written  as  o^/(l)  . In  the  theory  of  laser 

3-9 

speckle  patterns  , the  ratio  = o^/ (I  ) is  treated  as  a measure  of 
the  contrast  of  the  spatially  random  pattern.  Cj  = 1 for  Gaussian 
light.  It  is  convenient  to  extend  the  notion  of  the  contrast  of  I to  that 
of  the  energy  W received  by  the  detector.  The  contrast  of  W, 

C-yy  = aw/(W  ),  generalizes  and  can  be  interpreted  as  the  contrast 

of  the  measured  fading. 

The  general  formula  for  the  signal-to-noise  ratio,  rewritten  now 


yv-y-v  . 
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/ \l/* 

<nS  > 

SNR  = r , (8) 

f1+c2w<"s>+-5r^ 

yields  some  important  limiting  cases.  For  large  signals  (n  >>n  , C~z  ) 

o IN  W 

the  SNR  is  signal -fading -limited  and  equal  to  the  inverse  of  the  contrast: 

-l 

SNR  = C . For  noiseless  W,  (C  = 0),  the  SNR  depends  only  on 

w w 

signal-shot-noise  and  background  and  dark  current  noises: 

SNR  = (n  >/[(n  >+  (n  ) . For  background  noise  larger  than  the  signal, 

O O IN 

(n^  » n^),  the  SNR  is  background-noise -limited:  SNR  = (ng>  / (n^)  ^ . 

2 . 3 Fading -Limited  SNR 

This  section  discusses  the  fading  contrast,  C^.,  for  complex 
Gaussian  fields.  We  first  treat  fading  contrast  in  terms  of  the  temporal 
properties  of  the  illuminating  field  and  the  illuminator-object-imaging 
system  geometry,  then  apply  the  results  to  computations  of  fading  contrast 
for  optically  rough  and  smooth  surfaces  and  for  limiting  values  of  integration 
time  and  detector  aperture.  Several  practically  important  cases  of 
illumination  are  considered:  case  (a)  treats  cross  - spectrally  pure  light; 
case  (b)  addresses  monochromatic  illumination;  case  (c)  considers  the 
effects  of  polychromatic  illumination  on  the  contrast. 

To  gain  a qualitative  insight  into  computations  of  contrast  of  fading, 
consider  a laser -generated  field  at  the  object,  characterized  by 

some  coherence  area  and  time,  At  and  Tx.  The  field  leaving  the  object, 

V^\  has  coherence  area  and  time,  A2  and  T2.  If  the  target  is  stationary 


RIVERSIDE  RESEARCH  INSTITUTE 


with  respect  to  the  illuminating  beam,  the  coherence  times  of  V ^ and 


may  be  assumed  equal.  However,  the  transverse  coherence,  whose 


area  is  assumed  to  be  larger  than  the  object's  dimensions  , can  be 


modified  by  the  object  surface  roughness,  or  even  destroyed  completely. 


Propagation  of  the  field  V in  space  and  the  spatial  filtering  imposed 


by  the  imaging  telescope  lead  to  a further-modified  coherence  area  and 


time  of  the  image  field,  V . For  spatially  incoherent  objects,  for 


example,  the  coherence  area  A3  is  that  of  the  diffraction-limited 


resolution  cell  at  the  image  plane  ’ . The  coherence  time  T3  is  equal 


to  T2  if  the  maximum  differential  time  delay  in  propagating  along  different 


optical  paths  from  the  object  points  to  the  image  plane  is  much  smaller 


than  T2;  otherwise  T3  is  smaller  than  T2. 


2.3.1  General  Theor\ 


Computation  of  the  fading  contrast  is  a difficult  task  in  general; 


it  depends  on  the  covariance  function  of  the  image  irradiance  and  there- 


fore requires  evaluation  of  the  fourth  moment  of  the  field  at  the  image 


plane: 


C2  = 
W 


fdtJdtJd^jW  COV  { Ij  , I2} 


Jdtl  JdtZ  Id2xilld2xi2  <X1>  <X2  > 


The  integrations  in  time  and  space  are  over  T and  A respectively, 


Ij,  I2  are  irradiances  at  x.  , t3  and  x.  , t2,  and 


covf  ij  t i2}=  <ij  I2)-<I1)<I2>  . 


For  complex  Gaussian  fields!^’ ^ the  covariance  of  the  irradiance 


r ~ • * • • ■ i - - - 'v  - I,  —»*■»■ 
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can  be  expressed  in  terms  of  the  degree  of  coherence,  Y12(At). 

cov{l1(I2}=  1 Yiz(tll  t2)  |2  <I2  > <I2  >,  (H) 

where 

vi2(t1,t2)=  (V(3)(xil,t1)  V(3)>:!(x.2,t2))/((I1)(I2  ))lU.  (12) 

When  such  fields  are  stationary  and  linearly  polarized,  the  fading 
contrast  reduces  to 

CW  = A ”2  T "2  Jdtj  Jdt2  Jd2xu  fd2Xi2  | Y12(At)|2  (13) 

where  At  = |tj  - t2|.  The  reciprocal  of  this  is  often  interpreted  as 
the  number  of  degrees  of  freedom  in  linearly  polarized  light^. 

Using  Wolf's^  theory  of  propagation  of  the  coherence  function, 
the  degree  of  coherence  at  the  image  plane,  Yi2(At),  can  be  expressed^  in 
terms  of  the  degree  of  coherence  of  the  field  leaving  the  object,  Y^^At): 

Y12(At)=  F12  * < yj2  ^At  + R/c))  , (14) 

in  which  the  deterministic  transfer  operator  F12  is  applied  to 

(2 ) 

the  expected  value  of  the  propagated  Y12  (At).  R (x  Qi , x Q2  , x . xi2)/c  is 

the  differential  time  delay  in  propagating  along  different  optical  paths, 
one  from  xQl  to  and  another  from  xQ2  to  x^  . The  expectation  is 

taken  with  respect  to  realizations  of  R,  which  depends  on  random  variable 
z,  and  of  the  random  reflection  coefficient,  r]  . This  operation  is  dictated 
here  by  defining  the  coherence  to  apply  to  ensembles  of  both  light  sources 
and  scattering  surfaces  - 

If  the  differential  time  delay  R/c  is  much  smaller  than  the  coherence 
time  T2,  then  (vj2  \&t  + R/c))  can  be  factored  into 
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(Y^(At  + R/c)>  = (yj  ^ (At))(exp(j2TT v R/c)>  . (15) 

The  first  factor  is  the  expectation  of  the  degree  of  coherence  of  the 
emergent  field  at  the  object,  with  respect  to  the  random  reflection 
coefficient.  It  may  be  expressed  as  the  product  of  the  degree  of 
coherence  of  the  illuminating  field,  Y^(At),  and  the  mutual  correlation 
coefficient  of  r\  at  two  object  points,  P^Ax).  In  the  second  factor, 
the  path-length  difference  R can  be  decomposed  into  the  sum  of  a 
deterministic  part,  D(9),  and  a random  part,  Z(0).  The  former  depends 
on  the  angle  of  incidence  0 of  the  illumination  and  on  the  optical  path 
lengths  between  points  on  the  median  scattering  surface  and  the  image 
plane.  The  random  part, 

Z(9  ) = [z(x  Q1)  - z(x  02  )]  (1  + cos  0),  (16) 

involves  the  spatial  structure  of  the  surface  roughness,  zj^)  , and  the 
obliquity  factor. 

The  factor  due  to  the  surface  roughness  is  expressible  as 

(exp(j2rrv  Z/c)>  =0(k',  -k'),  (17) 

where  0(klt  k2)  = ( exp[j (k^  + k2z2) ])  is  the  joint  characteristic  function  of 
z (x  ) and  z(xn.)  and  k'  = 2tt(v/c)(1+  cos  0).  The  remaining  factor  combines 
with  the  deterministic  impulse  response  in  a convolution  that  accounts  for 
all  the  geometrical  effects,  including  the  telescope  aperture,  leaving 

V12(M)=.fd2xolJd2Xo!Hl!(*o..So2)«0''.  -k'>Pfl<45o>Y®W>.  (18) 

where  H12  is  the  impulse  response  in  the  absence  of  surface  roughness 
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and  for  uniform  reflection  properties, 
of  interest  can  now  be  examined. 


The  consequences  for  several  cases 


a.  Cross -Spectrally  Pure  Illumination 

If  in  addition  to  R/c  < < T2  = T,  the  spectral  density  of  the  light 
at  the  target  is  uniform,  the  last  factor  in  (18)  is  reducible  to  the  product 
of  the  spatial  and  temporal  degrees  of  coherence  y/^O)  and  Y^W)17. 
Thus,  for  Gaussian  fields  satisfying  the  condition  of  cross  spectral  purity, 
the  degree  of  coherence  at  the  image  plane  becomes  y (At)  = y (0)  y (At) 

12  12  u 

where 


V°>=Id2xoi  .rd2xo2Hi2(xoi>io2^(k'.  -k')PT1(Axo)Yi(21)(0)  (19) 

and  Yn(At)  = yn  (At).  This  indicates  that  the  degree  of  spatial  coherence 
at  the  image  plane  results  from  propagation  of  the  transverse  coherence 
of  the  illumination,  modified  by  the  object  surface  structure.  In  particular 
f°r  Yi2  ^ Vi2^°)  at  the  imaSe  depends  only  on  object  surface  structure 

and  on  the  coherent  point  spread  function  of  the  system.  On  the  other 
hand  the  temporal  degree  of  coherence  remains  unchanged  provided  that 


the  illuminator,  object,  and  receiver  do  not  move  with  respect  to  each 

other  during  the  measurement.  For  this  case  the  square  of  the  fading 

contrast,  C , is 
a 

Ca  =(T'ZJdt1  Jdt2|yi(i1)(At)|2)(A-2Jd2x.i  Jd^lvJ2)  (20) 

b.  Monochromatic  Illumination 

Laser  illumination  is  often  idealized  as  fully  monochromatic 
radiation: 


**,  •*,  r.  ^ •“  »*  .*  » _ «*  * * ,*  "_*•  * * “ • ' » " m " • r * ' » - » r » - » ~ 
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V teoi-^U&o)  exP[-j2rrVJ-  (21) 

It  follows  that  the  degree  of  coherence,  Y^(At),  is  given  by: 

12 

Y1(21)(^t)=  PTJ(Ax0  ) exp(-j2nvoAt),  (22) 

where  is  the  mutual  correlation  coefficient  of  the  amplitude  at 

two  object  points.  The  fading  contrast,  Cb,  is  then  time- independent: 


c^A‘2JdSi  J^-ia  lvi2l; 


c.  Polychromatic  Illumination 

Case  (b)  can  be  generalized  to  include  frequency  diversity  techniques 
which  decrease  the  contrast  of  the  speckle  pattern.  Let  the  target  be 
illuminated  with  M monochromatic  laser  lines,  of  equal  strength  but 
different  frequencies.  The  illumination  field  is 

m M 

VM  = u(-o)rJ=l  exP[-i2nvmt]  (24) 

(3) 

The  instantaneous  image  field,  , is  then  a superposition  of  M 


contributions, 


VmW> 


/ o \ 2 

Upon  integrating  the  instantaneous  irradiance  distribution,  I.  =|v'  M , 
over  a time,  T,  long  enough  to  average  out  all  beat  frequency  terms, 
the  irradiance  becomes  a sum  of  irradiances,  one  for  each  laser  line, 


and  the  energy  deposited  on  the  detector,  Wx , , is 
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Assuming  that  all  contributions  are  of  equal  strength  and  stationary, 


the  square  of  the  contrast,  Cw,  is 

M 


c2m=m_2a‘2K*i,  PS*  2 Iy  | 


m,  n=l  12 , mn 


where,  by  a simple  extension  of  (18), 


V = Jd2x  f d2x  H 
12. mn  J oi  J 02 


<»  Hl=,  mn  <°>  • (28| 


If  the  laser  lines  are  dense  enough  to  allow  replacement  of  the  summation 


over  lines  by  integration  over  frequency,  the  energy  W becomes 


Wc  = J d2xi  Jdv  S(v)  ^i-  V), 


where  S(v)  is  the  strength  of  illumination  at  frequency  V,  attenuated 


by  the  reflection  coefficient,  r]  , which  is  assumed  to  be  frequency- 


independent  and  uniform  across  the  object.  The  spatial  integrand  is  the 


image  intensity  distribution  from  quasimonochromatic  illumination  of  the 


object,  with  a coherence  time  that  does  not  satisfy  the  condition  for 


cross-spectral  purity  and  is  much  shorter  than  the  integration  time. 


This  model  for  the  image  intensity  has  been  used  to  study  polychromatic 


18  -21 

speckle  patterns.  Postulating,  for  simplicity,  uniform  S(v)  between 


Vi  and  v2  , the  square  of  the  contrast  is 


Cc  = B"2A’2J'd2xi2  J*  dvm J*  2<ivn  Gi2;i 

Vi  Vj 


where  B = V2-Vi  is  the  overall  bandwidth  and  G12,m  is  the  Fourier 


transform  of  yl2  mn17  . 


F*  - • • . ’ •*  fc  1 , k%V  . <0  I*  •*  • ^ m ► ft  » * t)\ 
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2.302  SNR  for  Rough  and  Specular  Surfaces 

The  main  results  of  the  previous  section,  Eqs.  (20,  23,  27,  30),  are 

valid  when  the  field  at  the  image  plane  is  a spatially  complex  Gaussian 

22 

process.  It  has  been  3hown,  experimentally  by  Ohtsubo  and  Asakura 

23 

and  theoretically  by  Goodman  , that  this  assumption  is  valid  for 

sufficiently  rough  surfaces.  Considerable  effort  has  been  devoted 

recently  to  a study  of  the  effects  of  surface  roughness  on  the  contrast 

of  monochromatic  ’ ^ and  polychromatic  speckle  patterns^ 

Qualitatively,  for  smoother  surfaces  the  random  phase  modulation  of  the 

field  leaving  the  object  is  distributed  over  a range  narrower  than  ± TT, 

the  corresponding  field  at  the  image  can  not  be  treated  as  a complex 

Gaussian  process,  and  the  contrast  of  the  speckle  pattern,  measured  as 

a^/  ( I ),  is  smaller  than  unity,  decreasing  as  the  surface  roughness 

decreases  . At  the  limit  of  specular  surfaces  the  joint  characteristic 

function  of  surface  roughness  approaches  unity.  For  monochromatic  and 

polychromatic  illumination  the  corresponding  energies  deposited  at  the 

detector  become  deterministic  and  the  contrasts  are  zero.  In  case  (a) 

the  contrast  depends  only  on  the  temporal  behavior  of  the  field. 

Surface  roughness  is  treated  here  as  a stationary,  zero-mean 

Gaussian  process.  The  joint  characteristic  function  required  for 

computation  of  y and  y is  given  by 

12  I2;mn 


0(ki,  k2)=  exp[-  ~ cr*  <k*  +k|  + 2ktkz  p)], 


(31) 


where  p(ixQ)  is  the  normalized  correlation  of  the  surface  profile 
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P(Axq)=  <z(x01)z(x02)>/a^  . (32) 

To  gain  more  qualitative  insight  into  the  effects  of  surface  structure 
on  the  degree  of  coherence  of  the  image  field,  let  the  correlation 
function  p(Ax<J)  fall  quadratically  from  unity  near  the  origin,  with 

*■ 

width  a;  p«  1 - | Ax^J  /a2.  Then  the  required  characteristic  function's 

h 

dependence  on  Ax0  is,  approximately, 


0(k',-k')~  exp[-(k'o  )2|  Ax  | /a2]. 


Takai  has  shown  that  the  above  approximation  for  a Gaussian  correlation 

is  accurate  within  a few  percent  for  k'oz£  TT . Furthermore,  for  large 

k'a  , the  function  can  be  approximated  by 
z 


0(k',-k)«TT  (-^-) 


and,  similarly, 


4(km-  -^)=“P[-T^,km-kk)2'k'  k'  ? 


m n z 


Physically,  these  approximations  imply  that,  for  rough  surfaces  (oz  2 X ), 

the  structure  behaves  as  if  it  were  composed  of  point  scatterers,  and 

25 

the  object  is  spatially  incoherent.  This  model  was  used  by  Goodman 
2-7 

and  many  others  to  study  speckle  pattern  statistics. 

Let  *W  and  °W  denote  the  energies  deposited  on  the  detector  for 
Images  of  rough  and  smooth  surfaces,  respectively.  The  corresponding 
contrasts  for  cases  a,  b and  c will  be  discussed  for  different  ratios  of 


1 - - - *■  - j ' . . 1 

• ■ W - V 'J.  - * • - • 
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coherence  time  and  coherence  area  to  the  integration  time  and  detector 
aperture.  Consider  first  *W  under  condition  (a)  of  cross  spectral 
purity.  Substituting  (34)  for  0(k',-k')  in  (19),  we  obtain  the  well-known 
Van  Cittert-Zernike  theorem  ^ that  the  spatial  degree  of  coherence 

i; 

equals  the  autocorrelation  function  of  the  impulse  response: 

*v„<°>  ' <*„,*„).  (36) 

as  appropriate  for  an  ideal  diffuser.  With  the  spatial  and  temporal 
coherence  separated,  the  fading  contrast  becomes 


*<■ 


1 


Al  Jj_ 
A T 


for  A3  > > A,  T3  >>  T 
for  A3  <<  A,  T3«T 


(37) 


where  Mandel's  definitions  of  coherence  area  A3  and  coherence  time 
T3  are  employed: 


a3  = Jd2x  |v12(x)|‘ 

T3  = Jd t I Yn(t)  I2 


(38) 

(39) 


For  specular  surfaces,  0(k',-k')=l  and  the  degree  of  spatial 
coherence  is  given  by 

M) 


Vl2  “ Id  X01  Jd  XQ2  HI2^2Eoi  >i5o2^Pri  ^5o  ^ ^12 


(40) 


For  uniform  illumination  (Y^  =1)  and  uniform  reflection  (P^  =1),  it 


follows  that 


12 


- 1.  Then  the  fading  contrast  is  given  by 


% \v.  y - V V Y '.-  V V-'a  .-.-. 
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°c 2 = <> 

a s 


1 


li 
^ T 


(41) 


lor  Tj>>  T 
for  T3<<  T 

For  specular  surfaces,  fading  becomes  minimal  when  the  coherence  time 
is  very  short  with  respect  to  the  integration  time,  regardless  of  the 
detector  aperture. 

b.  In  the  case  of  monochromatic  illumination,  the  degree  of  coher 


ence  for  a rough  scatterer  is  only  spatial  and  the  contrast  is  given  by 


*c£  = i 


i 

A, /A 


for 


A3  > > A 


A,  <<  A 


(42) 


For  specular  surfaces  °W^  is  deterministic  and 


°cb  = 0 


(43) 


c.  Lastly,  we  return  to  the  case  of  polychromatic  illumination. 

The  general  expression  (27)  for  can  be  evaluated  for  rough  surfaces 

(3) 

by  substituting  (35)  into  (19)  for  Y12.mn.  Assuming  that  the  wavelength 

27 

dependence  of  the  PSF  is  negligibly  small  , and  that  adjacent  lines  are 
separated  in  frequency  by  Av,  the  contrast  is  expressible  as 


*cl  _ , M-l 


# 2 _ 


M 


[*1  + 77  Z)  (M-n)exp(-bzn2)  T 
L M n=l  J 


(44) 


where  b = 2tt  (^v  oz/c)(l+ cos  0 ).  This  predicts  that  the  contrasts  for  the 
limiting  cases  of  complete  correlation  ( all  lines  have  the  same 
frequency)  and  large  decorrelation  (large  Av)  are  equal  to  and  to 

'''Cjij/M,  respectively.  The  result  can  be  simplified  for  the  following 


limited  cases: 
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CM  “ < 


1 Av  = 0 and  A3  > > A 

1/M  Av>c/az  and  A3  » A 
(A3/A)/M  AV>c/az  and  A >>  A3 


For  the  specular,  deterministic  case,  °C2  = 0. 

M 

For  polychromatic  illumination  of  continuous  spectrum,  the 
contrast  is  obtained  from  (30),  using  (35),  yielding 


*^2 


C2  = 
c 


“2  f dv  f dv  exp[-(k'  -k'  ) cr2  ] 
y J m J n ' m n z J 

vi 


(46) 


or 


21 


*C*e  = *cb['Tlet»-  ] (47) 

where  3 = 2tt(B  rr  /c)(  1+  cos  9).  The  contrast  decreases  with  increasing 
bandwidth  B.  Again,  for  the  specular  case,  °C2  = 0. 

The  ratio  A3/A  in  the  limiting  expressions  for  the  contrast  reveals 
that  it  is  possible  to  improve  fading -limited  SNR  by  spatial  averaging 
over  many  coherence  areas;  the  method  will  be  discussed  in  more  detail 
in  Section  2.5.  For  cross- spectrally  pure  light  the  improvement  in 
SNR  can  be  obtained  by  temporal  averaging  over  many  coherence  times 
of  the  illumination.  For  polychromatic  illumination  a gain  in  SNR  can  be 
obtained  from  an  increase  in  frequency  separation  between  the  lines,  or 
an  increase  in  bandwidth  of  the  illumination.  As  an  example,  consider 
the  situation  presented  in  Fig.  2,  showing  three  distinct  kinds  of  image 
resolution  cells  observable  in  laser -illuminated  systems.  Cell  1 is 
specular;  cell  2 is  optically  rough  but  nearly  normal  to  the  line  of  sight; 
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cell  3 is  also  rough  but  lies  skewed  to  the  line -of- sight  and  has  a 
considerable  depth  R^ , much  larger  than  the  roughness  scale. 

Assuming  the  same  average  energy  from  all  three  resolution  cells, 
the  fading -limited  SNRs  evaluated  for  each  cell  are  different.  Cell  1 
does  not  fade  at  all  in  coherent  illumination,  but  image  cells  2 and  3 
will  fade  considerably.  The  cells  will  also  behave  differently  as  the 
coherence  time  of  the  illumination  is  shortened,  or,  equivalently,  as  the 
bandwidth  is  broadened.  Cell  1 will  still  be  free  of  fading ; the  fading 
of  cell  2 will  diminish  slightly;  the  fading  of  the  "deep"  resolution  cell 
3,  will  diminish  most  rapidly  because  the  controlling  dimension  of  this 
cell  is  not  the  roughness  but  the  depth  of  the  cell. 

Figure  3 presents  fading -limited  SNRs  for  discrete  frequencies 

(solid  line)  and  for  a continuous  spectrum  (dashed  lines)  vs  the  total 

illumination  bandwidth,  for  different  scales  of  roughness.  For  discrete 

multiline  illumination,  the  smallest  and  largest  frequency  differences 

between  lines  are  AV  = B/(M-1)  and  B respectively.  The  SNR  ranges 

from  1 to  M1^,  the  latter  being  obtained  when  the  sum  in  (44)  is 

negligible  compared  to  unity.  For  the  case  of  two  lines,  and  to  reduce 

this  sum  to  less  than  0.01,  the  bandwidth  B=Av  must  exceed  2.14(c/4TTa  ) 

z 

Thus,  to  decorrelate  two  speckle  patterns,  each  due  to  a single  line 
illuminating  a surface  with  az  = 1 Urn,  the  frequency  offset  cannot  be 
smaller  than  50  THz.  For  wider  bandwidths,  the  SNR  is  greater 

when  more  lines  are  used,  but  for  smaller  bandwidths,  the  SNR  is  seen 
to  be  larger  for  M=2  than  for  M=ll.  This  can  be  explained  by 
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examination  of  the  sum  in  (44),  which  has  M-l=10  terms,  only  one  of 

them  smaller  than  0.01.  There  are  10  terms  offset  in  frequency  by 

B/10,  9 terms  by  B/9,  8 terms  by  B/8,  etc.  It  follows  that  if  the 

total  bandwidth  is  ten  times  broader  than  the  minimum  Av,  a surface 

roughness  with  a = lum  mutually  decorrelates  all  eleven  speckle 

patterns  and  the  corresponding  SNR  is  equal  to  ,/ll  . Generally,  for 

a fixed  total  bandwidth  B and  surface  roughness  a , the  SNR  is  larger 

z 

for  a smaller  number  of  lines  than  for  more  lines  when  the  bandwidth 
is  less  than  the  one  that  decorrelates  the  speckle  patterns  corresponding 
to  the  smaller  number  of  lines. 

For  a given  degree  of  roughness,  the  SNR  for  a larger  number 

of  lines  is  well  approximated  by  that  for  a continuous  spectrum  when 

the  bandwidth  is  smaller  than  (M-l)c/2rra  , for  which  not  all  the 

z 

speckle  patterns  are  decorrelated.  For  wider  bandwidths,  the  SNR 
can  be  well  approximated  by  M1^2  and  (82/tt)^4  for  discrete  and  continuous 
spectra,  respectively.  The  significant  difference  between  the  two  wideband 
cases  is  that  the  discrete-spectrum  SNR  is  independent  of  bandwidth, 
while  the  continuous -spectrum  SNR  is  proportional  to  the  square  root 
of  the  bandwidth. 

Typical  material  roughness  may  vary  from  fractions  of  a 
micrometer  to  several  tens  of  micrometers.  For  a surface  with  az=lum 
illuminated  with  radiation  covering  the  visible  region  0.4  um-0.  7 um 
(B=321  THz)  the  SNR~  3.  An  Ar  laser  generating  4 lines  at 
.514,  .494,  .484,  and  .474  um  , was  used  in  the  first  experiment 
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testing  the  usefulness  of  the  frequency  diversity  technique  for  decreasing 
the  speckle  contrast^*.  The  experiment  has  shown  that  the  contrast  of 
a superposition  of  speckle  patterns,  each  due  to  a single  line  diffracted 
by  a rough  surface  with  az  > 10  um  and  each  with  unity  contrast,  is 
noticeably  smaller  then  unity.  Fig.  3 confirms  that,  even  for  the 
smallest  frequency  offset  between  the  lines,  B = 12.6  THz 
there  is  some  decorrelation. 

2 . 4 SNR  for  Low-Level  Signal 

Consider  first  the  case  when  there  is  no  fading  (C  = 0)  and  the 

w 

SNR  reduces  to  (ng  )/ [ (ng  ) + (n^  ) ] ^.  By  the  preceding  discussion,  this 
formula  is  valid  when  temporally  coherent  and  partially  coherent  illumination 
fields  are  reflected  from  specular  surfaces.  Fading  limits  the  SNR  only 
for  strong  signals,  since  it  adds  C^.  <ng  f to  the  noise  term  above;  see  (8). 
Figure  4 is  a typical  plot  of  SNR  as  a function  of  signal  (nc  ) for 
particular  levels  of  noise  (nN>  and  fading  contrast  C w . For  spectrally 
pure  and  for  monochromatic  illumination,  Czw  is  the  inverse  of  the 
number  of  degrees  of  freedom  in  the  field  across  the  detector.  This 
number  becomes  the  effective  number  of  spatial  modes,  A/A3,  or  the 
number  of  temporal  modes,  T/T3,  whenever  either  of  these  is  large, 
or  the  number  of  spatio-temporal  modes,  (A/A3)(T/T3),  when  both  are 
large.  With  no  background  light  or  fading,  the  SNR  is  signal-shot-noise 
limited  and  equal  to  the  square  root  of  the  number  of  signal  photo- 
electrons; this  is  the  middle  asymptote  in  the  figure  . The  background 
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noise -limited  asymptote,  for  weak  signals,  is  proportional  to  (nc  ) . 

The  asymptotic  fading -limited  SNR,  for  strong  signals,  is  independent 
of  (ng>  and  equal  to  the  inverse  of  the  fading  contrast.  The  graph 
shows  the  combined  effects  of  signal  photocounts,  background  and  dark 
current  counts,  and  fading  due  to  partial  coherence  in  determining  the 
detector  SNR. 

2 . 5 Trade-Off  of  Spatial  Resolution  and  SNR 

To  maximize  the  spatial  resolution  of  the  stationary  imaging  system, 
the  detector  aperture  is  normally  made  no  larger  than  the  diffraction- 
limited  resolution  cell  of  the  telescope.  In  many  situations,  it  may  be 
desirable  to  improve  the  SNR  at  the  expense  of  resolution.  Larger 
apertures  can  increase  the  SNR  by  reducing  fading  contrast.  The  nature 
of  this  trade-off  may  be  appreciated  by  considering  a uniformly  illuminated 
and  uniformly  reflecting  target  viewed  against  a uniform  background  and 
taking  as  the  reference  SNR  the  one  associated  with  the  output  of  a 
detector  having  an  aperture  that  matches  the  diffraction -limited  resolution 
cell  of  the  telescope. 

In  the  general  expression  (8)  for  SNR,  both  (n  > and  <n  T>,  the 

s N 

expected  numbers  of  signal  and  noise  counts,  are  proportional  to  the 

area  A of  the  detector  aperture.  On  the  other  hand,  C^.  is  inversely. 

2 8 

proportional  to  A,  whenever  it  can  be  approximated  by  A3/A.  Dainty 
has  shown  that  this  is  valid  when  A/A3  exceeds  5 and  that  C2  is  only 
slightly  larger  than  A3/A  when  1 < A/A3  < 5.  Furthermore,  for  spatially 
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incoherent,  well-resolved  objects,  the  diffraction-limited  resolution  cell3,  4 
and  the  coherence  area  at  the  image  plane,  A3,  are  the  same  size.  It 
follows  that  both  C^.  (n^)  and  (n^)/(ng)  are  approximately  independent 
of  the  detector  area,  so  that  the  SNR  in  (8)  varies  as  A1/2  . One 
concludes  that  an  improvement  in  SNR  can  be  traded  against  a loss  of 
resolution  proportional  to  the  square  of  the  enhancement  in  SNR,  or  else 
in  direct  proportion  to  it  when  only  linear  resolution  is  considered. 

2 . 6 Conclusions 

For  images  of  actively  illuminated  targets  detected  by  a photon 
counting  array,  the  SNR  is  affected  by  the  coherence  of  the  radiation, 
the  roughness  of  the  target  and  the  PSF  of  the  optical  system,  as  well 
as  by  the  aperture  and  integration  time  of  the  detector.  All  these  effects 
can  be  accounted  for  in  a compact  expression,  (8),  for  the  SNR,  by 
introducing  the  fading  of  the  signal  across  the  detector  and  the  fading 
contrast  C,^,  which  generalizes  the  speckle  contrast.  The  expression 
indicates  that,  for  increasing  signal  levels,  the  SNR  is  limited  first  by 
background  noise,  then  by  signal  shot  noise,  and  ultimately  by  the  fading. 

The  fading  contrast  involves  stochastic,  spatial,  and  temporal 
averaging  of  the  second  moment  of  the  irradiance,  or  the  degree  of 
coherence,  which  propagates  from  the  rough  target  to  the  image  plane. 

The  roughness  enters  by  way  of  the  joint  characteristic  function  of  the 
surface  statistics.  For  a resolution  cell  skewed  with  respect  to  the 
direction  of  observation,  it  is  the  depth  of  that  cell  that  plays  the  role 
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of  the  roughness  scale. 

The  expressions  obtained  for  the  fading  contrast  reveal  quantitatively 
the  efficiency  of  various  means  for  improving  fading -limited  SNR,  by 
temporal  averaging,  by  spatial  averaging,  or  by  increased  bandwidth  of 
illumination.  Plots  of  the  SNR  against  bandwidth  for  various  roughness 
scales  and  against  signal  strength  for  various  noise  and  fading  levels 
have  been  presented  and  it  has  been  shown  that  the  SNR  can  be  improved 
at  the  expense  of  spatial  resolution. 
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2.8  Figure  Captions 

Fig.  1 Geometry  of  the  imaging  system.  0 = object, 

T = imaging  lens,  I = image,  A = array  of  detectors. 

Fig.  2 Imaging  of  rough  target,  showing  various 
object  resolution  cells:  1 = specular,  2 = optically 

rough  but  nearly  normal  to  the  line-of-sight , 3 = optically 
rough  but  skewed  to  the  line-of-sight. 

Fig.  3 Fading-limited  SNR  for  discrete  frequencies 
(solid  lines)  and  for  a continuous  spectrum  (dashed  lines) 
vs.  the  total  illumination  bandwidth,  for  different  scales  of 
roughness,  assumed  greater  than  the  illumination  wavelength. 


Fig.  4 Typical  plot  of  SNR  vs.  signal  level,  <n  >, 

5 

with  its  asymptotic  approximations.  The  steepest  asymptote 

is  for  background-noise-limited  operation:  SNR  = 

1/2 

< ns>/<nN>  ' . The  middle  asymptote  is  for  signal-shot- 

1/2 

noise-limited  operation:  SNR  = <n  > ' . The  horizontal 

s 

asymptote  represents  target-fading  limitation:  SNR  = Cw  ^ . 
Noise  level  < nN>  = 9,  fading  contrast  Cw  = 0.1. 
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III.  MAXIMUM  ANGULAR  ACCURACY  OF  PULSED  LASER  RADAR 
IN  PHOTOCOUNTING  LIMIT 


To  estimate  the  angular  position  of  targets  with  pulsed 
laser  radars,  their  images  may  be  sensed  with  a four-quadrant 
.noncoherent  detector  and  the  image  photocounting  distribution 
processed  to  obtain  the  angular  estimates.  The  choice  of  esti- 


bv- 


mator  depends  on  the  statistics  of  the  photocounts,  which  re- 
flect those  of  the  energy  deposited  on  the  detectors  by  the  la- 
ser field  returned  from  the  target.  For  an  optically  rough  tar- 
get, this  is  a spatio-temporal  complex  Gaussian  process,  produc- 
ing photoelectrons  that  obey  a negative  binomial  distribution 
with  M degrees  of  freedom.  The  limits  imposed  on  the  accuracy 
of  angular  estimation  by  signal  and  background  radiation  shot 
noise,  dark  current  noise  and  target  cross-section  fluctuations 
are  calculated.  Maximum  likelihood  estimates  of  angular  posi- 
tions are  derived  for  optically  rough  and  specular  targets  and 
their  performances  compared  with  theoretical  lower  bounds.  These 
estimates  are  realizable  in  systems  with  automatic  gain  control, 
fast  enough  to  follow  the  random  fluctuations  of  the  target  re- 
turns. For  slower  systems  and  large  signals,  the  accuracy  of 
the  resulting  estimates  is  limited  by  target  fluctuations. 


Portions  of  these  results  were  presented  at  the  1975  meeting 
of  the  Opt.  Soc.  of  Am.  in  Boston,  J.  Opt.  Soc.  Am.,  65,  1204A, 
with  M.  King,  W.  Edelson  and  P.  Diament. 

Accepted  for  publication  in  Appl.  Optics,  with  P.  Diament, 
M.  King  and  W.  Edelson. 
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3.1  Introduction 

The  potential  of  considerably  higher  accuracy  than  that 
of  microwave  radars  motivates  interest  in  employing  laser  radar 
for  angle  tracking.  With  the  advent  of  high -power  lasers,  an- 
gle-tracking laser  radars  are  particularly  attractive  for  space- 
borne  applications  because  the  absence  of  atmospheric  turbulence 
allows  their  potential  for  high  accuracy  to  be  used  to  full  ad- 
vantage. In  certain  satellite  tracking  applications,  objects 
at  distances  of  megameters  must  be  positioned  to  decimeters. 
High-accuracy  tracking,  used  for  precisely  pointing  a communi- 
cation laser  beam  from  one  satellite  to  another,  promotes  effi- 
cient use  of  laser  energy  in  the  communication  link.  It  has 
also  been  pointed  out1  that  precise  knowledge  of  the  relative 
positions  of  two  spacecraft  performing  a rendezvous  maneuver  is 
often  impractical  to  obtain  from  ground-based  radars?  an  accu- 
rate laser  radar  located  on  board  the  spacecraft  could  serve  as 
a guidance  aid. 

In  an  angle-tracking  system,  a sequence  of  direct  measure- 
ments of  target  position  is  fed  into  a tracking  filter  which 
may  produce  refined  estimates  of  past,  present,  and  future  tar- 
get position.  Generally,  the  accuracy  of  the  refined  estimate 
is  proportional  to  the  accuracy  of  direct  measurements,  so  that 
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schemes  for  improving  measurement  accuracy  are  potentially  of 
great  significance. 

Maximum  likelihood  estimation  of  the  location  of  an  opti- 
cal object  has  been  studied  by  Saleh,2  who  considered  shot-noise 
limited  detection  of  the  images.  McGarty3  has  investigated  es- 
timation of  the  optical  position  of  the  centroid  of  the  image 
photocounting  distribution  (IPD).  The  properties  of  the  IPD 
have  been  analyzed  by  Amoss  and  Davidson4  in  the  context  of  de- 
tection of  weak  optical  images,  and  by  Helstrom, 5,6  who  consid- 
ered object  radiation  in  the  presence  of  uniform  background 
noise.  All  these  papers  have  considered  only  self-illuminated 
noncoherent  objects. 

This  paper  considers  the  estimation  of  angular  positions 
of  laser-illuminated  optically  rough  and  specular  targets  from 
their  images.  The  field  of  the  returned  laser  light  is  treated 
as  a spatio-temporal  stochastic  process  that  forms  the  image 
irradiance  distribution  (IID),  which  is  sensed  by  a four-quad- 
rant  detector  and  converted  into  the  IPD.  An  estimate  of  the 
target  angular  position  (TAP)  is  obtained  by  splitting  and  dif- 
ferencing the  IPD. 

The  accuracy  of  estimates  of  optical-image  TAP  based  on 
four-quadrant  detection  has  been  evaluated  in  the  past  by  adapt- 
ing conventional  microwave  radar  theory. 7-9  The  optical  an- 
alogue of  the  classical  microwave  beam-slitting  formula, 

a = g (A/D)/SNR  , (1.1) 
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indicates  that  the  RMS  error  a is  proportional  to  the  optical 
system  resolution  (A/d)  and  inversely  proportional  to  the  de- 
tected SNR,  through  some  function  g of  the  image  distribu- 
tion. 

However,  three  important  differences  suggest  that  the  accu- 
racy of  TAP  estimates  for  laser  radar  systems  may  not  be  the  same  as  for 
microwave  radar  systems  . First,  the  laser  system  can  employ  non- 
coherent (or  direct)  detection  of  the  target  return  to  extract 
the  angular  position.  Second,  the  quantum  noise  may  be  quite 
apparent  in  the  laser  system,  though  insignificant  in  the  micro- 
wave  system.  Finally,  the  laser  system  can  resolve  targets 
with  modest-size  telescope  apertures  while  the  microwave  system 
deals  primarily  with  unresolved  objects. 

This  paper  considers  targets  considerably  smaller  than  a 
diffraction-limited  resolution  cell  of  the  imaging  telescope. 

The  case  of  a target  that  extends  over  many  such  resolution 
cells  has  also  been  considered10  and  will  be  discussed  in  a 
forthcoming  paper.  The  TAP  is  estimated  directly  from  the  out- 
puts of  the  photodetectors,  which  sense  a single-look  IID. 

During  this  single  look,  the  target  does  not  translate  with 
respect  to  the  optical  line -of -sight;  target  rotation  is  per- 
mitted during  a look,  however,  so  long  as  the  noncoherent  re- 
flectance distribution  of  the  target  remains  unchanged,  within 
the  resolution  limit  of  the  imaging  telescope.  The  mean  and 
standard  deviation  of  the  single-look  TAP  estimates  are  derived. 
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and  the  limitations  imposed  by  laser  speckles,  photocounting, 
background,  and  dark  current  noise  on  the  accuracy  of  the  es- 
timates are  investigated. 

1 Mathematical  Idealizations 


In  order  to  concentrate  on  essential  considerations,  sev- 
eral assumptions  have  been  made  in  this  analysis.  First,  let 
the  uniformly  reflective,  optically  rough  target  be  illuminated 
by  a linearly  polarized  laser  beam  of  wavelength  A.  The  rough 
surface  height  profile  is  regarded  as  a spatially  stationary 
random  process  with  standard  deviation  ag  and  coherence  area 
pg  and  does  not  alter  the  polarization  of  the  incident  light. 
The  target  size  is  taken  to  be  smaller  than  the  coherence  vol- 
ume associated  with  the  laser  illumim tion. 

The  aberration-free  imaging  system  has  a loss-free  square 
pupil  of  side  D,  an  effective  focal  length  F,  a space-in- 
variant  coherent  point  spread  function  (PSF)  h,  and  is  sit- 
uated in  the  far  field  of  the  target  at  a distance  R from  it. 
The  areas  of  the  diffraction-limited  resolution  cells  at  the 
target  and  image  planes  are,  respectively,  AQ  = (AR/d)2  and 
A^  = (AF/D)2.  The  area  Aq  is  assumed  to  be  much  larger  than 
the  target,  whose  surface  contains  many  surface-coherence  areas. 


I-  fc  - % • * • K - . • » » - - „ • '»  • > * k * W * * * W * - \W  ' * . 

.-•v  ;.y;..y.yy  y y. y ^y-y-y-v 


RIVERSIDE  RESEARCH  INSTITUTE 


The  detectors  are  assumed  to  be  ideal  photon-counting  de- 
vices with  an  ideal  time  response,  identical  quantum  efficien- 
cies r|,  and  uniform  spatial  responsivity. 

„ 2 Noise  Source 

The  sources  of  noise  considered  in  this  analysis  are  tar- 
get fading,  signal  and  background  light  shot  noise,  and  dark 
current.  These  are  discussed  briefly  in  the  following  para- 
graphs . 

According  to  the  scalar  theory  of  diffraction,  the  spatial 
coherence  of  the  illumination  is  destroyed  by  surface  roughness 
when  as^  A. 11,1 2 In  that  case,  the  field  returned  from  the 
target  can  be  treated  as  a superposition  of  temporally  coherent 
fields  originating  at  spatially  independent  scatterers  on  the 
target  surface,  with  each  one  introducing  into  the  outgoing  wave 
a random,  uniformly  distributed  phase  shift.  Since  the  total 
number  of  scatterers  is  large  (approximately  equal  to  the  ratio  of  the 
illuminated  target  area  to  the  surface  coherence  area  p 13 ), 
the  central  limit  theorem  implies  that  the  field  returned  from 
the  target  is  a stationary,  zero-mean,  complex  Gaussian  process 
over  the  receiver  aperture  plane.11  In  general,  the  coherence 
area  of  this  field  can  be  obtained  from  tht  Van  Cittert-Zernike 
theorem.14  For  instance,  if  the  target  is  a sphere  of  radius 
r,  then  the  coherence  area  is  given  by  A2  a tt  . (0.6  AR/2r)2. 

It  follows  that,  under  the  assumptions  made,  the  receiver  aper- 
ture is  smaller  than  A_ . 
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When  the  target  rotates,  the  receiver  pupil  plane  field 
becomes  a spatio-temporal  complex  Gaussian  process1 7 with  co- 
herence time  T2  equal  to  the  time  it  takes  for  the  coherence 
area  A2  to  travel  across  a point  aperture.  Therefore  the 
coherence  times  of  the  pupil  plane  field  and  of  the  correspond- 
ing image  field  are  equal. 

It  follows  that  the  instantaneous  image  irradiance  distri- 
bution i(x^jt)  is  a spatio-temporal  process,  obeying  a neg- 
ative exponential  probability  density  function  (PDF)  with  pa- 
rameter <(  l(x^,t)^>  , the  mean  irradiance. 

The  energy  deposited  on  a detector  of  area  A during 
time  T is : 


M=/dX  f 

A 7T 


dt  I ( x i , t ) 


(1.2) 


In  this  analysis  we  refer  to  the  fluctuations  in  W as 
the  target-fading  noise.  The  deposited  energy  W is  assumed 
to  obey  a gamma  PDF  with  parameter  y and  M degrees  of  free- 
dom. If  T3  is  the  coherence  time  of  the  image -plane  field 


A3  is  its  area. 

the  number 

of 

degrees  of 

freedom 

is  : 

M = (T/T3)(A/A3) 

for 

T » 

T3 

and 

A» 

a3 

(1.2a) 

M = A/A3 

for 

T « 

T3 

and 

A» 

a3 

(1.2b) 

M = T/T3 

for 

T » 

T3 

and 

A« 

a3 

(1.2c) 

M = 1 

for 

T « 

T3 

and 

A« 

a3 

( 1 . 2d ) 
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Mandel  has  shown1 9 that  Eq.  (1.2a)  applies  well  for  partially 
coherent  Gaussian  light  that  satisfies  the  conditions  for  cross- 
spectral  purity.  Equation  (1.2c)  is  satisfied  by  the  noncoher- 
ent background  light.  Equation  (l„2d)  is  well  satisfied  by 
monochromatic,  linearly  polarized  laser  speckle  patterns, 1 2)1 3’ 20 
and  Eq.  (1.2b)  is  a good  approximation  for  spatially  integrated 
laser  speckle  patterns . 1 2 > 1 3 > 20 

Consider  also  the  return  from  a target  illuminated  with  a 
multiline  laser.  For  optically  rough  targets,  all  components 
of  the  fields  are  stationary  complex  Gaussian  processes  over 
space.  The  components,  taken  to  be  of  equal  strength  for  simplicity, 
can  be  treated  as  statistically  independent  if  the  frequency  off- 
sets between  the  lines  are  sufficiently  large  to  allow  the  rough 
surface  to  produce  spatially  decorrelated  speckle  patterns .x 2 > 22 > 2 3 
For  an  integration  time  long  enough  to  average  out  he  beat 
frequency  cross  terms,  the  resulting  energy  obeys  a gamma  PDF 
with  the  number  of  degrees  of  freedom  equal  to  the  number  of 
laser  lines. 

In  this  analysis  we  will  often  use  as  the  fading-limited  SNR 
the  ratio  of  the  mean  to  the  standard  deviation  of  W.  It  can  be 
shown  that  in  many  situations  of  interest,  this  ratio  is  equal 
to  the  inverse  of  the  square  root  of  the  number  of  degrees  of 
freedom. 12»13*18 

Spatially  incoherent  background  noise  with  a coherence 
time  much  shorter  than  the  integration  time  T is  assumed. 
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Photocounting  noise  results  from  the  random  fluctuations 
of  the  total  number  of  photocounts,  m,  at  the  detector's  out- 
put,, The  nuirber  m is  a sum  of  three  terms:  the  counts  gen- 

erated by  the  noisy  signal  W,  the  counts  due  to  the  back- 
ground radiation,  and  the  dark  current  counts.  The  three  terms 
are  mutually  independent;  the  first  one  is  a doubly  stochastic 
pocess,18  the  other  two  obey  Poisson  PDFs.  The  impact  of  the 
first  term  on  the  system  performance  is  of  particular  inter- 
est here. 

This  analysis  does  not  take  into  account  the  target-angle 
noise  that  may  result  from  random  fluctuations  due  to  finite 
size.  This  effect  is  caused  by  variations  in  the  phase  front 
of  radiation  from  a multiple-point  target  as  the  target  changes 
its  aspect.  It  has  been  found  that  in  microwave  monopulse  the 
effects  of  this  noise  decrease  with  the  target  range. 26 ’ 27 
Thus  the  angle  noise  can  be  disregarded  for  ranges  long  enough 
for  phase-front  distortions  to  be  negligible.  It  is  also  as- 
sumed that  detector  and  preamplifier  noise  have  been  reduced 
to  the  point  where  kT«  hv,  so  that  the  Johnson  noise  is 
negligibly  small.  Also,  1/f  noise  is  not  considered. 
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3 . 2 Problem  Formulation 

The  problem  of  TAP  estimation  is  addressed  systematically 
by  (1)  finding  the  theoretical  limits  given  by  the  Cramer-Rao 
lower  bounds  for  the  standard  deviations  of  the  estimates  and 
(2)  by  proposing  operations  on  the  data  that  lead  to  estimates 
with  standard  deviations  that  approach  the  theoretical  lower 
bounds.  In  this  paper,  the  Maximum  Likelihood  Estimates  (MLE) 
of  the  TAP  are  considered  as  the  candidates  likely  to  possess 
the  desired  properties . In  order  to  derive  the  Cramer-Rao 
bounds  and  the  MLEs,  knowledge  of  the  joint  probability  den- 
sity function  (PDF)  of  the  stochastic  processes  is  required. 

In  this  section,  the  PDFs  for  the  cases  treated  in  the  paper 
are  introduced.  They  are  used  in  Parts  3.3  and  3.4  for  the  de- 
rivations of  the  Cramer-Rao  bounds  and  of  the  MLEs. 

For  the  unresolved  target,  the  image  shape  matches  the 
noncoherent  PSF  of  the  imaging  system  and  the  location  of  the 
peak  of  the  I ID  depends  on  the  target  location.  Its  peak  value, 
and  therefore  the  total  energy  in  the  image  field,  can  be  treated 
as  a random  variable  that  depends  on  the  realizations  of  the 
field  in  front  of  the  receiver  aperture. 

The  analysis  considers  targets  that  are  displaced  by  an 
unknown  angle  in  the  x direction  with  respect  to  the  optical 
axis  of  the  imaging  system;  see  Fig.  1.  For  simplicity,  only 
one  of  the  two  angular  coordinates  is  treated  since,  for  a 
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square  aperture,  the  two  are  independent.  Also,  for  unresolved 
targets  the  fluctuations  do  not  depend  on  the  segmentation.  The  im- 
age irradiance  is  redistributed  between  the  two  photodetectors  by  the 
reflective  prism,  the  proportions  depending  on  the  TAP.  If  the  irra- 
diance is  a temporally  stationary  stochastic  process,  the  energies 
deposited  on  the  two  detectors  during  an  integration  time  T are 


Wx  = SXW  , W2  = s2w 


(2.1) 


where  is  the  energy  captured  by  the  pupil  aperture 

during  time  T.  The  fractions  Sx  and  S2  depend  on  the  noncoherent 

PSF,  i.e.,  the  squared  magnitude  of  the  coherent  h(xQ  ,yQ  ;x^,y^ ),  and 
are  given  by 


(2.2) 


It  is  assumed  that  the  PSF  is  normalized  to  fix  the  total 
energy  in  the  intensity  PSF  at  unity: 

Sx  + S2  = 1 • (2.3a) 

In  addition,  if  0O  is  smaller  than  the  diffraction-limited 
resolution  of  the  aperture,  then  approximately, 

si  ~ s2  * le0  • (2.3b) 

where  l is  a known  proportionality  coefficient  depending  on 
the  PSF.  For  the  square  aperture  under  consideration,  £ = 2D /X. 
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Let  the  number  of  photoelectrons  resulting  from  the  pho- 
tons arriving  at  the  two  detectors  in  random  quantities 
l^WSj/hv  = NSX  and  T}WS2/hv  = NS2,  be  nx  and  n2,  where  r) 
is  the  quantum  efficiency  of  the  detector,  and  hv  is  the 
energy  of  a single  photon.  N = r\\j/hv  obeys  the  gamma  PDF 

with  M degrees  of  freedom20  1 28  and  parameter  7 = (M/N  ): 

s 

M ~ 1 

P(N)  = 7 "(M  - ~l)l  eXp  ( -7N)  (2,4) 

with  N ^ 0 only.  The  limiting  forms  are: 

p(N)  = exp(  -N/Ng)/Ns,  M = 1 (2.4a) 

p(N)  = 6 (N  - N ) , M ->  00  (2.4b) 

s 

where  N = T]  < W )>/hv  is  the  mean  signal  photocount  collected 
by  the  pupil.  Generally,  as  the  number  of  degrees  of  freedom 
M grows  large,  the  relative  fluctuations  in  W become  small? 
limiting  forms  for  evaluating  M were  discussed  in  Sec.  3.1.2. 

For  the  analysis  that  follows,  it  is  useful  to  obtain  a 
formula  for  the  joint  PDF  of  nt  and  n2  by  using  the  Mandel 
transformation1 8 of  the  random  light  distribution: 


P(n1,n2)  = [ dN  P(N)p(nl,n2|N)  , (2.5) 

where  p(nx,n2|N)  is  the  joint  conditional  PDF  of  nx  and  n2 , 
given  the  probable  total  number  of  counts,  N: 
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P(nx,n2  | N) 


2 

n 

k=l 


-v 

V" 


(2.6) 


The  form  of  the  Mandel  transformation  given  in  (2.5)  uses  a 
common  value  of  N,  rather  than  N^Nj,  because  in  the  model 
for  unresolved  targets  the  energies  deposited  on  the  two  detec- 
tors are  derived  on  the  average  from  the  same  coherence  area  of 
the  pupil  plane  field. 

If  the  laser  radiation  returned  from  the  target  is  accom- 
panied by  background  noise  radiation,  the  photocounts  at  the 
detector  outputs  are  the  sums  of  the  signal  and  background 
counts.  It  is  assumed  that  the  background  photocounts  at  each 
detector  obey  a Poisson  PDF  with  a mean  count  N^/2 . The 
detector  dark  current  noise  can  also  contribute  to  the  total 
number  of  counts  at  the  detector  output.  It  is  assumed  that 
the  dark  current  counts  of  each  detector  obey  a Poisson  PDF21 
with  an  average  count  N^/2 . Thus,  the  total  counts  in  each 
detector,  m1,m2,  are  the  sum  of  the  mutually  independent 
signal,  background,  and  dark  current  counts.  The  noisy-case 
joint  PDF  of  mx  and  m2  is 


q(ml,m2)=l  dN  p(N  )q  (n^  ,m2  |N  ) , 

J0 


(2.7 


•— 


$ 

v 


where  q(m1,m2|N)  is  the  joint  conditional  PDF  of  and  m2 
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q(m  ,m2  |N  ) = II 
k=l 


2 -(SkN+NN/2) 

n.  s — sn <SkN+V2> 


(2.8) 


where  NN  = (N^+N^).  Equation  (2.7)  can  be  obtained  directly 
from  (2,5),  for  the  sum  of  mutually  independent  Poisson  random 
variables  is  a Poisson  random  variable  with  mean  value  equal 
to  the  sum  of  the  means . 2 4 

The  formula  (2.7)  can  be  made  more  suitable  for  deriving 
the  MLE  of  the  TAP  if  it  is  approximated  by  the  expression  valid 
for  small  background  and  dark  current  noises  (N  < 1)  and  for 
large  signal-to-noise  ratio  ( Ng/NN » 1 ) . Using  in  this  for- 
mula the  approximations 

(SkN  + NN/2)mk  = (SkN)mk[l+  (i\/Sk)(NN/2N)] 


exp(  -Nn)  = 1 - Nn  , 

the  latter  serving  merely  to  maintain  the  proper  normalization, 
the  approximate  joint  PDF  can  be  written  as: 


m + m S mi  S,m2 


where 


+m2  -H 


< N e~w  > - 1 +mx  + m^ 

NM  m + m - 1 “ l M + N_  )NS  ' 

< N 2 e N N V S ' 


(2.9) 


(2.10) 


m ' •*«.'  • ''i  * ■ ^ . - . V"  . * . . . . 
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and  the  expectations  are  taken  with  respect  to  N,  using 
(2.4). 

For  the  special  case  when  the  random  fluctuations  of  N 


are  negligibly  small  (number  of  degrees  of  freedom  M is  large), 
the  PDF  of  N can  be  well  approximated  by  the  Dirac  delta  func- 
tion28 in  (2.4b).  After  letting  the  number  of  photocounts  for 
this  case  be  m and  m , the  joint  PDF  can  be  obtained 
from  (2.7)  and  (2.8)  as 

= l(,"0.>,n02ll,S)  • (2‘U) 

Note  that  the  above  joint  PDFs  are  functions  of  the  TAP,  6q , 
via  the  parameters  Sx  and  S2. 

The  first  problem,  considered  in  the  next  section,  can  be 
stated  as:  given  the  measurements  nx  and  n2  (for  a fluc- 

tuating target  in  the  absence  of  background  and  dark  current 
noise)  and  the  a priori  knowledge  of  Ng;  or,  given  the  mea- 
surements mQi  and  mQ  2 (nonfluctuating  target  in  the  pres- 
sence  of  background  noise)  and  the  a priori  knowledge  of  Ng 
and  NN,  what  are  the  corresponding  theoretical  lower  bounds 
for  the  standard  deviations  of  the  TAP  estimates? 

The  second  problem,  addressed  in  Part  3.4,  can  be  stated 
as:  given  nx , n2  and  Ng,  or  given  m1,m2  (fluctuating 

target  in  the  presence  of  background  noise)  and  Ng,  NN,  and 
Mr  or  given  m0i,  mQ2  and  Ng  and  N^,  (l)  find  the  expres- 

sions for  the  parameter  60  that  maximize  the  joint  PDFs 
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p(nx.»n2)j  or  q0  (m0lJm02),  respectively;  (2)  cal- 

culate the  means  and  the  standard  deviations  of  these  expres- 
sions. It  is  assumed  that  the  parameters  Ng,  NN,  and  M are 
known  as  a result  of  averaging  many  previous  measurements,  or 
otherwise . 


„ - W.  - *4  V „“>•  _ ■ - JN  - *-  L%L>  . ‘.'.M*  - > S*  * > .. 
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(3.3) 
are  given 

(3.4) 

< nxn2>  = sxS2n|(1  + £)  . (3.5) 

The  inequality  for  the  lower  bound  is  therefore 

var  9*  (l~2  -0j)/Ns  . (3.5) 

Since  GQ  has  been  assumed  to  be  much  smaller  than  l/l , the 
lower  bound  for  the  standard  deviation  of  the  estimate  becomes 
effectively 

°min  - VW.1/*  • (3.7) 

The  expression  for  the  lower  bound  of  the  standard  devia- 
tion of  9*  is  the  same  as  the  beam  splitting  formula  given  by 
(1.1 ),  with  the  SNR  equal  the  square  root  of  the  signal  count 
Ng  and  l~1  proportional  to  A/D.  This  implies  that  the  theo- 
retical limit  for  the  accuracy  of  the  TAP  derived  from  nx  and 
n2,  given  Ng,  is  limited  by  signal  shot  noise  and  does  not 
depend  on  target  fading. 


< 


Sin  p(n  ,n  ) )2 

> ■ <*v«>  < 


ft  • a 


> 


The  expected  values  / n£  )>  (k  = 1,2)  and  < nin2  /> 
by 

<"k>  = <Nssk>2<1  + S>  + s*Ns  ’ 
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3.3.2  Nonfluctuating  Targets  with  Noise 

For  nonfluctuating  targets  the  inequality  (3.1)  can  be 
rewritten  as: 

var<^><  [Slnqo(mol,mo2)/S0o  ]2  >-1  , (3. 

where  <f>*  is  an  estimate  of  0Q  and  qQ  (mQ  x , mQ  2 ) is  given 
by  (2.11)  and  (2.8). 

Following  the  same  steps  as  above,  one  has 


2 

(Sin  q (in  .m  ) \ An  m \ 

"Se~-0,2,  ) > = 


> A 


(3.9) 


where  Nk  = SkNg  + nn/2  . With 


Ojk>  - (Nk  + ijNj 


< m m „ > = N N 

x oi  o 2 / 12  3 

the  inequality  (3.8)  can  be  written  as 


nq  + nm  r 

var  K > 5|  - eo 


Ji-T 

ns  + hnJ 


Again  neglecting  0Q  compared  to  l/l , the  lower  bound  of  the 


standard  deviation  is 


'o  min 


■ (ns+V,/2/'w,s  • 


• _*•  - . . 
A A ■ . ■ 

v‘  *-  a' 
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This  is  the  same  as  the  beam-splitting  formula,  with  the  SNR 
now  given  by  the  ratio  N /(Nc+N  )1/2.  This  suggests  that  the 
theoretical  limit  for  the  accuracy  of  the  TAP  derived  from  mQ x 
and  mQ2,  given  Ng  and  N^,  depends  on  the  shot  noise  gen- 
erated by  the  target-returned  laser  illumination  and  the  back- 
ground radiation,  and  on  the  dark  current.  It  should  be  noted 
that  (3.13)  applies  also  as  a bound  for  the  accuracy  of  angular 
position  estimates  for  temporally  noncoherent  targets,  for 
which  M is  large. 


The  theoretical  bounds  for  the  standard  deviations  of  es- 


B 

f. 


r ’ * 


timates  of  the  TAP,  such  as  those  obtained  above,  do  not  iden- 
tify the  estimator  whose  standard  deviation  attains  these 
bounds.  The  choice  of  the  MLEs  as  good  candidates  for  estimat- 
ing the  TAPs  is  made  here  because  they  are  asymptotically  effi- 
cient estimates;  i.e.,  their  variances  attain  the  Cramer-Rao 
bounds  as  the  number  of  statistically  independent  measurements 
increases . 

3.4.1  MLE  FOR  FADING  TARGET  WITHOUT  BACKGROUND  (1^  M<  oo,  = 0 ) 

The  MLE  of  the  TAP  0Q , is  that  value  of  0O  that  is 
most  likely  to  cause  the  observed  values  of  nx  and  n2  to 
occur.  It  maximizes  the  likelihood  function  p(nx,n2|0o)  for 
given  nt,n2.  If  this  maximum  is  interior  to  the  range  of  0O 
and  the  function  is  appropriately  continuous,  the  MLE  is  found 
from  the  likelihood  equation25 


r: 
kv 
/ ' ■ 


dp(nt,n  I 0Q )/d0o  = 0 . 


(4.1) 


By  (2.5),  (2.6)  and  the  fact  that  0Q  appears  only  in  Sl,S2, 


3 


ES- 

3 
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n n _ n +n  _ . . 

0 0 


ni‘  n2 


nx ! n2 ! ^0 


/ nx  na\ 

: (s>  s> ) = 


and  d 


/ ni  n2\  /N 

lsi  S2  J/aeo  -*-s  solved  by 

S = l~1(ni  ~ n2  Mni  +n2)  . 


(4.2) 


The  processing  required  to  realize  the  MLE  of  TAP  requires 
taking  the  scaled  ratio  of  the  photocount  difference  to  the 
total  number  of  photocounts,  where  the  scaling  factor  l is 
known  and  depends  on  the  shape  of  the  PSF  of  the  imaging  tel- 
escope. This  estimator  is  conceptually  similar  to  the  one  em- 
ployed in  microwave  angle -tracking  monopulse  systems.  9 In  the 
event  that  both  nx  and  n2  are  zero,  it  is  intended  that  the 
processor  not  provide  any  TAP  estimate  at  all,  as  if  no  light 
were  incident. 
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The  derivation  of  the  mean  and  variance  of  the  MLE  of  the 
TAP,  0 , is  presented  in  the  appendix  but  the  main  results 


are 

are 

stated  and  discussed  here.  The  mean  and  variance  of 
given  by  (A. 9)  and  (A. 18): 

So 

<®0  > = e„(i-J'M) 

(4.3) 

and 

var  0Q 

= -t_2J_MG  + 02J_M(1  - j"M -G) 

(4.4) 

where 


J = 1 

G = In  J 

For  the  strongest  fading 
rewritten  in  a form  more 

<§o> 

var  S0  = l~2  (in  Np  )/fcp  - 92q  [ ( In  Np  )/fcp  - Ng/N*  ] (4.8) 

In  ( Nq  + 1) 

^ 2 5 tA  o — » 

where  N = N_  + 1 and  the 
P 


Ng  + 1 ^ 

last  form  accounts  for  0 <<(  l~x . 


+ 1/y  = 1 + N_/M 

O 

(4.5) 

M-l 

+ Z-f  (J™  - l)/m  . 

m=l 

(4.6) 

(M-l)  the  above  equations 

can  be 

useful  for  interpretation: 

- floV<NS  + 1> 

(4.7) 

A - *• 

'/// 

, * a > a *a  ’a. 
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For  small  signals  and  strong  fading  the  bias  cannot  be  disre- 

A 

garded.  Figure  2 plots  <\Q0'}/Q0  as  a function  of  Ng  for 
different  values  of  M,  For  M = 1 more  than  50  signal  photo- 
counts are  required  to  reduce  the  bias  to  the  same  extent  as 
can  just  four  signal  counts  for  M^>  5. 

The  variance  of  0Q  can  be  written  as  a sum  of  terms 
dependent  on  and  independent  of  the  target  position  0q.  Both 
terms  of  var  0Q  are  dependent  on  the  target  fading  and  both 
can  be  made  negligibly  small  for  sufficiently  large  Ng . Fig- 
ure 3 is  a plot  of  the  RMS  error29[var  0Q  ]^2,  in  A/D  units, 
as  a function  of  N for  different  values  of  M.  The  region 
covered  is  that  of  0 small  with  respect  to  the  diffraction 
angle  A/D.  A square  DXD  pupil  aperture  was  assumed  for  which 
'i  = 2D/A0  In  the  small  photocount  region  2 ^ Nc  10,  target 
fading  does  not  play  a dominant  role  and  the  accuracy  of  §0  is 

limited  by  the  signal  shot  noise.  For  larger  Nc  the  limita- 

b 


\ 


RIVERSIDE  RESEARCH  INSTITUTE 


tions  imposed  by  target  fading  are  clearly  pronounced.  An  in- 
crease in  the  number  of  degrees  of  freedom  M is  particularly 
advantageous  for  fixed  Ng  when  the  fading  is  strong  (M  small). 

Although  0Q  is  biased,  it  is  instructive  to  compare  its 

standard  deviation  with  the  lower  bound  a . derived  for 

min 

unbiased  estimates.  The  lower  bound  is  plotted  on  Fig.  3 as  a 
dashed  line.  From  this  comparison  one  concludes  that  for  weak 
fading  and  for  strong  signals,  for  which  the  bias  becomes  neg- 
ligible, amin  can  be  used  as  a good  approximation  for  determ- 

A 

ining  the  accuracy  of  0 . Figure  3 indicates  that  this  approx- 
imation would  be  very  poor  for  strong  fading.  Suppose  that  an  RMS 
error  no  larger  than  0.1  \/D  is  desired.  To  achieve  this  accuracy 
the  expected  number  of  signal  photocounts  should  be  larger  than 
120  for  M = 1 and  larger  than  25  for  M)>  10.  Therefore,  to  obtain 
an  accuracy  of  0.1  k/D,  the  required  laser  power  must  be  five 
times  larger  for  M = 1 than  for  M)>  10. 

In  those  instances  when  the  bias  is  not  negligible,  the  ques- 
tion of  whether  to  modify  the  estimator  to  remove  the  bias  can 
readily  be  treated  by  comparing  the  accuracies  of  both  types  of 
estimates,  using  the  results  presented  here. 

Several  consequences  for  the  design  of  illumination  wave- 
fronts  and  measurement  schemes  result  from  the  above  analysis  of 
the  variance  of  the  MLE.  Qualitatively,  the  target  fading  lim- 
itations can  be  decreased  by  decreasing  the  temporal  coherence 
of  illumination,  as  illustrated  by  the  following  three  cases. 
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1.  A stationary  rough  target  is  illuminated  with  a sin- 
gle pulse  from  a multiline  laser.  The  illumination  wavefront 
is  a superposition  of  M laser  lines  mutually  offset  in  fre- 
quency. The  frequency  offsets  are  assumed  to  be  sufficiently 
large  so  that  the  speckle  patterns  at  the  receiver,  each  due 
to  a single  line,  are  mutually  decorrelated. 

2.  A spinning  target  is  illuminated  with  M short, 
equal -strength  pulses  from  a single-line  laser.  The  return 
from  each  pulse,  rendered  statistically  independent  from  the 
others  by  the  spin,  is  used  to  obtain  a MLE  of  the  TAP.  All 
M estimates  are  used  to  obtain  the  sample  mean  of  the  esti- 
mates of  the  TAP.  The  RMS  error  of  the  sample  mean  is  smaller 
than  that  of  a single-pulse  estimate  by  a factor  of  M1/2  . 

3.  A spinning  target  is  illuminated  with  a long  pulse 
from  a single-line  laser.  During  the  integration  time, which 
is  not  shorter  than  the  pulse  duration,  let  an  average  of  M 
speckles  be  captured  by  the  receiver  aperture. 

Cases  1 and  3 are  formally  similar:  a single  measure- 

ment is  made,  based  on  the  photocounts  generated  by  energy  W 
with  M degrees  of  freedom.  In  Case  2,  M measurements  are 
taken,  each  based  on  the  photocounts  generated  by  W with  one 
degree  of  freedom.  Combining  Cases  1 and  3,  by  use  of  multi- 
line  radiation,  one  can  perform  a single  measurement  in  which 
W is  sensed  with  an  effective  degree  of  freedom  equal  to  the 
product  of  those  of  the  two  cases.  Combining  Cases  1 and  2, 
by  use  of  M pulses  of  M lines,  M measurements  are  per- 
formed, each  resulting  from  sensing  W with  M degrees 


, " • ' 
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of  freedom.  All  three  cases  and  their  combinations  provide 
diversities  instrumental  in  easing  target-fading  limitations. 

Given  that  the  total  signal  Ng  is  fixed  in  all  three 
cases,  is  a smaller  RMS  error  obtained  from  a single  measure- 
ment with  M degrees  of  freedom  and  Ng/M  photocounts  per 

degree  of  freedom  (Cases  1 and  3 ) or  from  M measurements  each 
based  on  Ng/M  total  photocounts  with  one  degree  of  freedom 

(Case  2)?  Generally,  the  RMS  error  of  the  MLE  obtained  from 
a single  measurement,  based  on  energy  with  M degrees  of  free- 
dom, is  smaller  than  the  one  obtained  from  M independent  mea- 
surements, each  based  on  energy  with  one  degree  of  freedom,  pro- 
vided that  other  conditions  are  unchanged.  In  particular,  this 
implies  that  measurements  based  on  the  long  illumination  pulse 
offer  better  results  than  the  one  based  on  averaging  over  many 
short  pulses. 

Two  numerical  examples  are  presented  as  an  illustration  of 
the  use  of  the  owt„(N_,m)  formula,  based  on  (4.4). 

MLi-h  b 

Example  1 ; Let  the  total  signal  be  Ng  = 100.  Dividing  the 
signal  equally  among  10  laser  lines  that  are  sufficiently  off- 
set in  frequency  to  produce  10  degrees  of  freedom  in  the  energy 
deposited  on  the  detector,  we  obtain  an  RMS  error  (100,10) 

= 0.053  A/D.  If,  instead,  the  total  signal  Ng  = 100  is  divided 
equally  among  10  pulses,  the  RMS  error  of  a single  measurement 
is  c„__(10,l)  = 0.23  A/D  and  the  RMS  error  obtained  by  averag- 
ing  10  measurements  is  0.074  A/D. 
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Example  2 : For  the  same  total  signal  Nc  = 100  and  dividing 

the  signal  equally  among  5 laser  lines,  one  obtains  for  the 
RMS  error  aMLE(100,5)  = 0.056  A/D . The  RMS  error  of  the  MLE 
obtained  by  averaging  over  5 pulses  each  with  signal  Ng  = 20 
is  aMLE( 20,l)/5^2=  0.085  A/D. 

In  some  situations,  it  may  be  simpler  to  replace  the  de- 
nominator of  Eq.  (4.2)  by  <(nx  + n/> , an  estimate  of  this  ex- 
pectation being  obtained  by  averaging  the  total  number  of  counts 
over  many  observations,  each  representing  an  independent  reali- 
zation of  the  process.  The  resulting  estimator  is  not  a MLE. 

It  is  expressed  by: 

^ n - n 

0i  = 4-i  —1 L-  . (4.10  ) 

<ni  +n2> 

Formally,  the  expectation  in  the  denominator  makes  it  statisti- 
cally independent  of  the  numerator.  A similar  situation  occurs 
in  microwave  amplitude  monopulse  radars  with  slow  AGC:  the  time 

period  over  which  the  denominator  is  averaged,  the  integration 
time,  is  long  with  respect  to  the  inverse  of  the  target  scintil- 
lation bandwidth,  so  that  the  sum  channel's  fluctuations  are  elim 
inated  while  the  difference  channel  remains  stochastic.  It  has  been 
found  that  the  accuracy  of  tracking  of  microwave  systems  is  limited  by 
target  scintillations  (target  fading)  for  slow  AGC  and  free  of  these 
limitations  for  fast  AGC . 2 6 * 2 7 our  results  demonstrate  that  the  same 

A A 

is  true  for  optical  systems,  with  and  0Q  analogous  to  the  esti- 
mators obtained  in  microwave  monopulse  radar  with  slow  and  fast 
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AGCs,  respectively.,  The  performance  of  0X , compared  with 

A 

that  of  Qq , is  discussed  below. 

The  mean  and  variance  of  estimator  are 


<el> 

= ('e,d)-1<n1  -n2> 

(4.1 

(td)-2 

[°n. +0L  -2  )j 

i 

(4.1 

where  d = <^nx  + n2)>,  cov(n1>n2)  is  the  covariance  of  n^ 
and  n2,  and  and  c22  are  the  variances  of  nx  and  n 

The  photocounts  nx  and  n2  obey  the  marginal  PDFs  p(nx), 
P(n2)>  which  can  be  obtained  from  the  joint  one  in  Eq.  (2.5) 


as 


dNp (N ) 


(S^N)"* 

nk! 


(4.1: 


The  mean  and  the  variance  of 

<"k> 


are  given  by29 
= SkNS 


(4.1^ 


°nk  = SkNS  + SkNl/M  (4-^ 

and  the  covariance  is 

cov(nx,n2)  = S1S2N|/M  . (4.1€ 
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Substituting  these  into  (4.11,  4.12),  yields  the  mean  and 
variance  of  61 : 


Thus,  0X  is  an  unbiased  estimator  of  0Q  . Its  variance  is 
a sum  of  two  terms;  the  first  depends  on  the  PSF  of  the  system 
(via  t ) and  is  inversely  proportional  to  the  square  of  the  sig- 
nal-shot-noise-limited SNR.  The  second  term  is  proportional 
to  9*,  and  is  inversely  proportional  to  the  square  of  the 
fading-limited  SNR.  For  a target  on  the  optical  axis,  dQ  =0, 
the  fading  noise  can  be  disregarded  and  the  system  performance 


<0i>=0o  - 


var  0X  = 'E'“2/Ns  + 0J/M 


(4.17) 


(4.18) 


is  limited  by  the  signal-shot-noise.  For  0Q  ^ 0 and  large 

N the  photocounting  noise  can  be  neglected  with  respect  to 
S 

the  target-fading  noise  and  the  second  term  limits  the  accuracy 

of  this  estimator.  For  a fixed  and  nonzero  value  of  0Q , this 

term  is  largest  for  the  deepest  fading  (M  = l).  As  M increases, 

the  second  term  decreases  and  can  be  neglected  at  the  limit  of 

noncoherent  light  (large  M). 

Comparing  the  RMS  error  of  0Q  and  0X  with  the  lower 

bound  a . for  the  case  of  strong  target-fading  noise  (M  = l, 
min  * 

0Q=O.l  A/D),  one  concludes  that  (l)  the  theoretical  limit 
for  the  RMS  errors  does  not  depend  on  the  target  fading,  whereas 

A A f 

the  RMS  erros  of  0Q  and  0,  are  target  noise  dependent; 
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(2)  the  RMS  error  of  6 , which  is  an  unbiased  estimate,  is 
larger  than  for  aH  values  of  signal;  they  are  equal 

only  if  9q  = 0 or  if  M is  very  large.  The  difference  be- 
tween them  is  negligible  for  small  signals.  For  large  signals 
the  RMS  error  of  0 is  significantly  larger  than  a -n,  being 
three  times  larger  than  a . for  N„  = 10 3 and  a hundred 
times  larger  for  Ng  = 104?  (3)  the  estimate  9Q  is  biased, 

which  allows  its  RMS  error  to  be  smaller  than  a . for  small 

min 

A 

signals.  For  Ns  = l,  for  example,  for  which  the  bias  of  9Q 

is  50  per  cent,  the  RMS  error  of  0 and  o . are  0.3  a/d 

o min 

and  0.5  A/D,  respectively.  The  RMS  error  of  is  larger 

A A 

than  that  of  9X  for  values  of  Ng  for  which  0Q  is  approx- 
imately unbiased  and  signal  shot  noise  dominates.  As  the  sig- 

A 

nal  increases  to  where  the  target-fading  noise  dominates,  9Q 

• • A A 

performs  significantly  better  than  0X : the  RMS  errors  of  9Q 

are  2.4  and  6.6  times  smaller  than  the  corresponding  values  for 
0 for  Nc  = 10 3 and  Nc  = 104,  respectively. 

loo 
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discussed,  based  on  computer-simulated  data.  Confidence  in 
the  simulations  is  gained  from  their  good  agreement  with  the 
analytically  derived  formulas  for  M = 1 and  M = <». 

The  MLE  is  the  value  of  0Q  that  maximizes  qi(ml,m2) 
given  in  (2.9).  Differentiation  yields  the  condition 


VS1  SJ(  N>  [ s*  s|  J2n„ 


(4.19) 


But  to  zero  order  in  the  noise  N, 


N J 


mx  m 

S7  ~ S7  “ mx  + m2 

Applying  these  equations  to  the  factor  of  NN  in  (4.19),  the 
maximum  likelihood  condition  becomes,  to  first  order  in  N^, 

mxS2  -m2s1  = (mx  - m2  )NN/2NM  . (4.19a) 

Combining  this  with  (2.3a)  and  (2.3b),  one  finds  that  the  noise- 
less MLE  is  modified  by  the  factor  1 + N../N. . : 

N M 


A 


o 


l-1 


m - in 

1 , - 1 
mx  + m2 


M + NS 

M - 1 + mx  + m2 


(4.20) 


To  realize  A it  is  necessary  to  know  N„,  N„  and  M. 
These  parameters  are  assumed  to  be  known  from  averaging  of  many 
independent  measurements,  of  the  background  alone  to  obtain  an 
estimate  of  N^,  of  the  target-returned  laser  illumination  with 
the  background  to  obtain  estimates  of  the  sum  Ng+NN,  and 
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r.  »•. 
y.  ■ . 


u» 


lastly  of  the  variance  of  the  total  number  of  photocounts, 
var(mi+m2).  All  these  estimates  are  needed  to  estimate  M 
as 


M 


= N|[var(mi  +m2  ) - (Ng  + N^)  ] " , 


-x 


(4. 


since  var(mx  + m2  ) = Ng  + NN  + var  N and  var  N = N|/M. 


To  calculate  the  mean  and  variance  of  A it  is  conven- 

o 


ient  to  introduce  the  following  notation. 


Ao  AS  + AN  * 


= <(qs  + qn)  e_N> 


(4. 

(4. 


where  subscripts  S and  N are  used  to  indicate  terms  that 
are  independent  and  dependent  on  the  background  noise,  respec 
tively,  and  the  expectation  is  with  respect  to  N;  i„e.. 


as  = AsK,ma)  = i'1  s-rnjr  - 


X 2 


AN  = NN(1+MAS)AS("'1^2)/(M  -l  + mx  +m2)  , 


rox  , .m2 


fSxN)  (S2N) 

QsK*"a)  = — sn 


m2  ! 


°N  = NN^QS^mi,m2  " 1)+iQg(mi  " l,m2)  -Qs(mx^m2^ 


The  approximate  expressions  for  the  mean  and  variance  of  A , 


to  first  order  in  NN/Ng,  are 


<A„>  = <§„>  + <e'N^(QsAN+QNAs)>  , 


(4, 


.-■.••.'V-V. 


RIVERSIDE  RESEARCH  INSTITUTE 


75 


var  Aq  = var 


K + < e'NZVs <As  ‘ 2 < So  > » + 2<  (As‘<S0>N> 


(4.25) 

where  the  sums  are  over  and  m2 , exclusive  of  m1=m2=0J 

the  expectations  are  taken  over  N.  and  -C'  0 \ and  var  0 

x o / o 

were  computed  in  Sec.  (3.4.1).  The  remaining  terms  in  (4.24-25) 
can  be  evaluated  by  extension  of  the  computational  techniques 
employed  in  the  Appendix.  The  expressions  involve  modified  ex- 
ponential integrals  and  are  not  informative  in  general.  For 
M = 1,  however,  the  analytic  expression  for  the  mean  of  the 
estimator. 


<0  = 0 (NQ  + Nj/(NC  + 1) 


(4.26) 


valid  for  smaller  than  both  and  unity,  shows  that  the 

noise,  when  weak,  reduces  the  bias  of  the  MLE,  because  the  noise 
and  signal  bias  act  in  opposite  directions. 


• v v »-  .*  v > r • v v\- v v / v v * - 

- - • , ■*  «■«'..*  . ■ - _ * . ■ B •»  _ -r.  _ « • _ » • • « , »•  . • _ ■ •»  - _■ 


. * • • r * - - • “ * ' , 

. s. v.  \r/, 
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Standard  deviations  of  the  MLE  A were  obtained  in  com- 

.0 

puter  simulations  using  5000  realizations  of  the  random  mx 
and  m2  and  examined  as  functions  of  the  signal  Ng  for  dif- 
ferent values  of  M and  . For  large  NgJ  the  RMS  error 

A A 

of  Ao  can  be  considered  practically  equal  to  that  of  9q . 

For  weak  signals,  however,  the  difference  between  the  errors 
of  these  estimators  can  be  appreciable.  The  percentage  dif- 
ference between  the  RMS  errors  of  A and  6 , normalized  to 

o o 

A 

that  of  9 , is  plotted  in  Fig.  4 as  a function  of  the  signal 
Ng,  for  M = 1,2,20,  and  for  = 0.1  (solid  lines)  and 
Nn  = 0.5  (dashed  lines).  The  RMS  error  of  Aq  is  larger  than 
that  of  9q  for  the  same  N^,  and  M0  The  normalized  er- 

ror difference  decreases  with  the  signal  Nc  and,  for  fixed 
Ns  and  N^,  decreases  with  M.  Also,  for  fixed  M and  Ng, 
it  increases  with  the  noise  level. 

Consider  now  an  estimator  Ax  realizable  with  a slow  AGC 
system.  This  estimator  is  defined  as: 


(4.2' 


The  mean  and  variance  of  Ai  are  readily  obtained  from  those 
of  mi ± m2  from  (2.8),  without  approximation. 


var  A = 
x 


= fVNs  + xV/<Ns  + V ■ 

K . ,-2  HS  + »,1  p»S  + 

M NH  N_  + N.. 


(4.28 

(4.2S 


as 


The  likelihood  function  for  this  case  is  qQ  (mQ  x ,m0  2 | Ng  ) , 
in  (2.7a),  and  the  MLE  is  most  easily  obtained  by  maximizing 
In  qQ „ This  estimator  is 


A 


= I'1 


(4.30) 


The  mean  and  variance  of  $Q  can  be  obtained  by  the  method 
used  in  the  Appendix. 


a„>  - e„[1-fe'v]  > (4-31) 

-N  f -N  "1  -N 

var  0O  = eje  T|l-e  T + Ein  (-NT ) - l~2  (N^/Ns  ) 2e  TEin(-NT)  , 

(4.32) 

where  N = Nc  + N.T  and  the  exponential  integral  is  the  one  that 
appears  in  (A-13).  The  MLE  is  biased,  but  underestimates  0o 
only  slightly  whenever  the  total  of  signal  and  noise  is  not 

A 

negligible.  In  this  case  the  RMS  error  in  0Q  can  be  consid- 
ered equal  to  the  lower  bound  oQ  m^n  f°r  all  practical  purposes 
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Finally,  consider  an  estimate  $ , 


*1 


Nt 

£-1  _1 


in  - m „ 
0 1 0 2 


N_  <m„  , + m ) 
S x 0 1 02/ 


From  the  mean  and  variance  of  m ± m , 

01  02’ 


<*!>  - »0 


var  =<.-2(Ns+NN)/N| 

Thus,  this  estimator  is  unbiased  and  its 
equals  the  lower  bound  a 

o mm 


defined  as 

(4. 

there  follows  that 

(4. 

(4, 

standard  deviation 
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SUMMARY 

In  a laser  radar  monopulse  system,  the  angular  position 
of  a target  is  estimated  by  measuring  the  location  of  the  tar- 
get's photocounting  image  by  means  of  a noncoherent  four-quad- 
rant  detector.  The  limitations  imposed  by  target  fading  and 
by  signal  and  background  shot  and  dark  current  noise  on  the  ac- 
curacy of  the  estimate  have  been  investigated.  The  theoret- 
ical limits  given  by  the  Cramer-Rao  lower  bounds  for  the  stan- 
dard deviations  of  the  estimates  have  been  derived  for  fading 
and  nonfading  targets.  For  several  cases  it  has  been  shown 
that  the  lower  bounds  are  given  by  the  beam-splitting  formula; 
i.e.,  that  the  standard  deviations  are  proportional  to  the  dif- 
fraction-limited angular  resolution  of  the  imaging  telescope 
and  inversely  proportional  to  the  SNR,  where  the  noise  includes 
the  shot  noise  generated  by  the  target  and  background  radiation 
and  dark  current,  and  do  not  depend  on  the  target  fading. 

Maximum  likelihood  estimates  have  been  derived  for  fading 
and  nonfading  targets  and  their  means  and  standard  deviations 
have  been  evaluated.  All  these  estimates  share  the  following 
properties : 

1.  They  depend  on  the  familiar  monopulse  ratio  which, 
for  the  photocounting  problem,  is  given  by  the  ratio  of  the 


VL^rvrv  v xn  kt 
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difference  of  the  total  number  of  photocounts  at  the  detector 
outputs  to  their  sum.  This  normalization  makes  the  angular 
measurement  independent  of  the  target  range  and  significantly 
reduces  the  impact  of  target  fading  on  its  accuracy. 

2.  They  are  biased.  For  nonf luctuating  targets,  the 
bias  decreases  exponentially  with  the  signal  and  background 
noise.  For  fading  targets  the  bias  depends  on  the  number  of 
degrees  of  freedom,  M,  and  decreases  monotonically  with  the 
signal  at  a rate  that  is  significantly  faster  for  M y 5 than 
for  M = 1 . For  M y 5 and  N 2,  the  estimates  for  fluctu- 

O 

ating  targets  can  be  considered  practically  unbiased. 

3.  The  standard  deviations  of  the  estimates  approach  the 
theoretical  Cramer-Rao  limits  for  unbiased  estimates  when  the 
signals  are  large  enough  to  make  the  MLEs  unbiased. 

4.  The  processing  of  data  according  to  the  MLEs  requires 
systems  with  fast  automatic  gain  control,  which  can  follow  the 
temporal  fluctuations  of  the  target  fading.  The  fast  AGC  is 
necessary  to  obtain  the  correlation  between  the  difference  and 
sum  channels.  If  the  AGC  is  slow,  the  difference  and  sum  pho- 
tocounts are  decorrelated  and  the  corresponding  estimates, 
although  unbiased,  have,  standard  deviations  that  are  limited 
for  large  signals  by  target  fading  noise. 
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APPENDIX:  MEAN  AND  VARIANCE  OF  MLE 


We  need  <r)>  and  <r2>  for 


r = 

n1+n2  ■ 


where  the  counts  nx,n2  are  doubly  stochastic: 


with 


00 

<'*>  -J  z rs  p(n1,n2  |N)p(N)dN 
0 nx'n2 


(S1 N)  1 (S2N)  2 

p(niJn2  In) ^ s-p-  » 


Ti^nr  --7N 


(A.l) 


(A. 2) 


(A.  3) 


(A. 4) 


The  sum  in  (A. 2)  must  exclude  the  possibility  of  no  counts  at 
all,  na  = n2  = 0.  Let, 


ni  n2 

n - n n +n„  S S 


En  -n  n +n,  S S 
-hr  - 2 x 1 2 -L_  -1 

ni+n2  nx!  n: 


(A. 5) 


ni'n2 


E(xs 
(nx  ~n2)  — — 


(xsj  1 (xs  ) 


, (A. 6) 


ni,n- 


,'w 

.\  . • ■ v"i  %•  - 

■ ■ . --  V.  -■  v „■ 


Then 
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which  readily  sums  to 

(S+Sjx 

x^H/dx  = x(sx  - S2  ) e = x(Sx 

since  Sl+S2=l.  Integrating  this,  with 
H(x,S1,S2)  = (S4  -S2)(eX  -1) 
It  follows  that,  setting  J = 1 + 1/y, 


h(o,si,s2)  =0 


<r> 


oo 

= f H(N,S1, 


S2)  e_Np(N)dN 


= (Si  -S2)  <1  -e_N> 


= (S1  - S ) (1  - J"M) 


by  use  of  (A04) 


Similarly,  to  obtain  the  variance,  evaluate 

, 2 


ni  _n2 


nx,n2  ' 


by  forming 


a , 3k>  V , ,,(xSJ  ‘ <xS=>  2 

ST(x  & Zj  ( 1 '"a’  ! 


n2  I 


ni  jn2 


which,  after  rewriting  (nx  -n2)2  as 

[nl(n1  -1) +nx  - 2nxn2 +n2(n2  -1) +n2],  readily  sums  to 


m 

m 
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x — (x  — )=(x2S^+xS1-  2x2S1S  + x2S2  + xs2  )eX 

dx  <$x 


= [x2(S1  -S2)2 +x]eX  . 


(A. 12) 


Integrating  twice  yields  K in  terms  of  the  exponential  integral: 


x 

Kfx.S^S,)  = [1-  (S,  -Sa)a]  Je_U.=  l du+  (ex_1 


)(S1  -s2)2 


(A. 13) 


= 4S1S2[-Ein(-x)  ] + (S1  -S2)2(eX-l) 


Thus , 


N 

<r2>  = jK(N,S1,S2)e~N; 


p(N)dN 


(A. 14) 


= 4S1S2<-e"N  Ein(-N)>  + (S1  - S2  )2< l-e"“ > . 


The  double  integral  for  the  expectation  involving  the  exponen- 
tial integral  can  be  evaluated  as  follows: 


N 

00  m 

<-e‘N  Ein(-N)>  = f dNp  (N  )e_N  J 


00  00 

•'n  •'v 


dNp  (N)e 


m 


= (— l-)m 

+ yJ 


= (r^)M  + + 


00 

- I dx  xm  1 |exp(-yx)  - exp(-[l+7]x)J 


M (l+±)“-l 


. (A. 15) 


v'-.-'v . v-v-’ 


A % •• 
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Hence,  with  J = 1 + 1/7, 

var  r = 4SXS2J 
or,  since  43^  = 1 - (sx  - S2  )2 


var  r 


M-l 

G = InJ  + 

m=l 


-M  / 

J G + (Sx  - S2 


+ (sx  -s2)2j-m(i 

and  letting 

( J™  “ 1 )/m  , 

)2J-M(1  - J_M  -G) 


J-M) 
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Fig.  1. 


Fig.  2. 
Fig.  3. 

Fig.  4. 


Figure  Captions 

Single-Axis  tracking  sensor.  L = laser,  0 = object, 
T = telescope,  Q = four -quadrant  prism,  R = receiver 
irradiance,  I = image  irradiance,  nx,n2  = output 
photocounts . The  areas  under  the  irradiance  curve 
on  the  two  sides  of  the  center  line  are  in  the  ratio 

Si:S2. 

A 

Normalized  mean  of  MLE  0Q  vs  mean  signal  count. 

RMS  error  of  MLE  Qq  vs  mean  signal  count.  The 
dashed  line  is  the  Cramer-Rao  lower  bound. 

A 

Normalized  difference  between  RMS  errors  of  and 

o 

8 ilOCXu,  -ca)/oa  vs  mean  signal  count  for  noise 

0 ii  u u 

counts  N.=0.1  (solid  lines)  and  N..  = 0.5  (dashed 


WHS 


Angular  tracking  of  targets  with  monopulse  laser  radar  can 
be  performed  by  tracking  the  centroid  of  the  image.  The  image 
irradiance  distribution  is  sensed  with  an  array  of  noncoherent 
photocounting  detectors  and  the  individual  counts  processed  to  pro- 
vide estimates  of  the  image  centroid  coordinates  and  of  its  size  and 
orientation.  These  estimates  are  discussed  for  optically  rough 
and  smooth  extended  targets,  viewed  against  a noncoherent  uni- 
form background.  General  expressions  for  the  mean  and  variance 
of  the  image  centroid  and  second  moment  are  obtained  as  a func- 
tion of  coherence  of  illumination,  target  extent,  number  of 
degrees  of  freedom  in  the  target's  image,  the  diffraction-lim- 
ited resolution  of  the  image  telescope,  and  the  integration 
time  of  the  detectors.  The  fundamental  limits  imposed  on  the 
accuracy  of  the  estimates  by  the  shot  noise  from  the  image  and 
background  irradiances  and  by  the  laser  speckle  noise  are 
discussed. 

Portions  of  these  results  were  presented  at  the  1975  meeting 
of  the  Opt.  Soc.  of  Amer.  in  Boston,  J.  Opt.  Soc.  Amer.,  61, 
1204A. 

Submitted  for  publication  in  Appl.  Optics  with  P.  Diament. 
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4.1  Introduction 

An  important  application  of  laser  radar  is  to  angular  track' 
ing,  in  which  a sequence  of  direct  measurements  of  target  posi- 
tion is  fed  into  a tracking  filter  to  produce  refined  estimates 
of  present  and  future  target  position.  The  theory  of  angular 
tracking  developed  for  microwave  radars 1)2  does  not  apply  fully 
to  laser  radars.  Unlike  microwave  radars,  laser  systems  can  re- 
solve targets  with  only  modest-size  telescope  apertures.  The 
quantum  noise  may  be  quite  apparent  in  the  optical  electromag- 
netic spectrum.  If  background  radiation  noise  is  not  a dom- 
inant factor,  the  laser  energy  returned  by  the  target  can  be 
reliably  sensed  with  noncoherent  detectors.  These  are  simpler 
than  heterodyne  sensors  and  are  particularly  attractive  for  use 
with  recently  developed  powerful  chemical  laser  sources  that 
have  poor  temporal  coherence. 

This  paper  discusses  estimation  of  angular  position,  size,  and 
orientation  of  extended  targets  by  use  of  monopulse  laser  radars 
which  derive  their  information  about  target  angular  position 
(TAP)  and  target  size  and  orientation  (TSO)  from  a single  pulse. 
The  TAP  is  estimated  from  the  position  of  the  centroid  of  the 
image.  The  image  centroid,  one  of  many  ad  hoc  estimators  that 
can  be  used  for  TAP  estimation,  generalizes  to  extended  targets 
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the  usual3  point  target4  angular  estimators  derived  from  four- 


quadrant  image  sensing  detectors.  The  TSO  is  estimated  from 


the  properly  normalized  image  second-moment  tensor. 


The  imaging  system  receives  both  reflected  light  from  the 


target^  which  is  illuminated  with  the  laser  beam  provided  for 


tracking,  and  from  uniform,  noncoherent  background  radiation. 


This  light  is  imaged  onto  the  detector  array,  which  converts 


the  photon  flow  into  an  image  photocounting  distribution  (IPD), 


the  electrical  signal  used  by  the  image  processor.  The  image 


processor  is  to  estimate  accurately  the  image  position,  size  and 


orientation  from  each  received  pulse. 


The  sensor  consists  of  an  array  of  L detectors;  the  k-th 


one  is  located  at  position  r^  and  registers  m^  photoelectron 


counts  in  response  to  one  laser  pulse.  Of  these,  some  are  back- 


ground photoelectrons  and  dark  current  electrons,  the  rest  are 


signal  photoelectrons.  In  the  absence  of  a target  image,  the 


noise  count  consists  of  only  background  and  dark  current,  fluc- 


tuating in  number  about  a mean  noise  count  < mN  )>  that  is  as- 


sumed to  be  known  and  the  same  for  all  the  detectors.  The  re- 


duced count  at  the  k-th  detector  is 


"k  = "k  - < mN  > 


the  total  count  diminished  by  the  known  mean  noise  count;  this 


includes  the  random  signal  counts  and  the  fluctuations  of  the 


noise  counts. 
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The  processor  considered  here  provides  the  tracking  system 
with  the  location  of  the  reduced  IPD  centroid  and  second  moment. 
The  former  is  the  vector  (see  Fig.  l) 

£ v2}‘  (2) 

and  the  latter  is  the  optical  moment  of  inertia  tensor 

£ ' (3) 

The  summations  are  over  all  the  detectors  of  the  array.  The 
moments  are  measured  with  respect  to  the  optical  axis  of  the 
imaging  system.  Upon  shifting  I to  refer  to  the  optical  cen- 
troid, its  eigenvalues  and  eigenvectors  provide  the  image  size 
and  orientation. 

The  numerators  of  these  quantities  estimate  the  first  and 
second  moments  of  the  image  irradiance  distribution  and,  except 
for  their  randomness,  correspond  to  those  of  the  irradiance  • 
distribution  at  the  target.  The  normalization  denominator  for 
both  moments  is  the  total  nuitber  of  counts  in  the  array,  reduced 
by  the  expected  value  of  the  total  number  of  noise  counts.  This 
normalization  factor  makes  the  position,  size,  and  orientation 
measurements  independent  of  the  target's  range  and,  because  it 
too  is  random,  can  ease  the  dependence  of  the  measured  moments 
on  the  target-related  random  fluctuations,  provided  a correla- 
tion exists  between  the  total  count  and  each  of  the  moments. 
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The  statistical  properties  of  the  estimators  depend  on 
those  of  the  IPD.  It  is  assumed  that  the  noise  counts  obey 
Poisson  statistics  with  parameter  equal  to  the  sum  of  the  mean 
background  and  dark  current  counts  and  are  independent  of  each 
other  and  of  the  signal  counts.  The  statistics  of  the  signal 
counts  reflect  those  of  the  energy  deposited  on  the  in- 

dividual detectors  by  the  laser  field  returned  from  the  target. 
These  counts  are  generally  correlated. 

As  has  been  shown,5  the  statistical  properties  of  the  en- 
ergy at  each  detector  depend  on  the  coherence  of  the  illuminat- 
ing light,  the  target-surface  roughness,  the  width  of  the  tel- 
escope point-spread  functions,  the  illuminator-target-telescope 
geometry,  and  the  aperture  and  integration  time  of  the  detec- 
tor. The  random  fluctuations  in  W make  the  resultant  signal 
counts  a doubly  stochastic  Poisson  process.6  For  an  optically 
rough  target,  the  returned  field  is  a stochastic  spatio-tempo- 
ral, zero-mean  complex  Gaussian  process  and  the  photoelectron 
counts  obey  a negative  binomial  distribution  with  M degrees 
of  freedom. 5)7  The  image  irradiance  distribution  from  such 
fields  has  the  characteristic  random  spatial  structure  referred 
to  as  the  laser  speckle  pattern. B> 9 For  an  optically  smooth 
target  and  temporally  coherent  illumination,  the  image  field 
does  not  fluctuate  and  the  photoelectrons  are  then  only  singly 
stochastic,  obeying  Poisson  statistics . 5> 7 
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To  determine  the  limitations  imposed  on  the  accuracies  of 
the  estimators  by  the  laser  speckle  structure,  the  signal  and 
background  shot  noise  and  the  dark  current,  the  general  ex- 
pressions for  the  means  and  variances  of  c and  I are  de- 
rived  in  the  next  section.  The  relevant  expressions  are  spe- 
cialized in  Sec.  4.3  for  TAP  estimation  with  the  monopulse 


laser. 
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MEAN  AND  VARIANCE  OF  THE 

The  performance  and  accuracy  of  the  estimators  q,  and  J 
are  determined  by  their  means  and  covariance  matrices.  These 
can  be  evaluated  by  adapting  the  method  of  statistical  error 
propagation10  or  the  delta  method11  to  the  centroid12  and  sec- 
ond moment  expressions.  It  is  shown  that  the  means  and  var- 
iances for  both  estimators  are  expressible  in  terms  of  a de- 
tection matrix  comprised  of  the  covariances  of  pairs  of  counts 

= <Cmk  " <mk>]K  ' <m-e>3>  • (4) 

Both  moments,  as  well  as  higl  er-order  ones  that  could  be 
treated  analogously,  are  of  the  form  of  a deterministic  func- 
tion of  the  location  of  the  k-th  detector,  averaged  over 

the  reduced  counts  at  each  detector: 

E ■ (5) 

This  random  quantity  is  intended  to  estimate  the  corresponding 
property  of  the  image 

F =2Z\<nk>/^<n<.>  • («) 

The  quantity  and  hence  also  E and  F,  can  be  of  any 

scalar,  vector,  or  tensor  character.  For  the  centroid  <J>  is 
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~kJ  E **■’  and  F •'■s  the  true  image  centroid  s.  For 

the  second  moment  tensor,  is  r^r^,  E is  I,  and  F 

is  the  true  optical  moment  of  inertia  tensor  M,  referred  to 
the  optical  axis . 

Let  N = £ be  the  total  mean  reduced  count  and  let 
Ax  = x - <^x)>  be  the  deviation  from  its  mean  of  any  random  var 
iable  x.  When  expressed  in  terms  of  such  deviations,  the  es- 
timator is 


E 


( F + ^\Ank/N )/  ( 1 +^"\nt/N ) 


(7) 


Expanding  the  reciprocal  of  the  denominator  in  a Taylor  series 
and  retaining  terms  to  only  the  second  order  in  An, 


e = f + ^2  (®k  ■ F^Ank/N  - ^2 


(\  - F)AnkAn^/^: 


(8) 


k,-t 


It  follows  that  the  mean  of  the  estimator  is,  to  this  order, 

< E > = F (4>k  -FXAr^An^/N2  (9) 

k,-t 


and  its  variance,  or  covariance  matrix,  is  given  to  the  same 
order  by  the  outer  product 


<AEAE> 


($k  - P ) (*£_  - F ) <AnkAn^ >/N2 


(10) 
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The  relevant  statistical  parameters  are  seen  to  be  only  the 
total  mean  N and  the  covariances  <(An^An^)>  of  the  reduced 
counts.  Note  that  the  estimator  is  generally  biased,  the  bias 
arising  from  the  correlation  between  the  random  numerator  and 
denominator  of  the  expression  for  E.  Note  also  that,  since 
the  difference  between  the  total  count  it^  and  the  reduced 
one  n^  is  not  random,  their  covariances  are  the  same: 
<AnkAnt>  = <AmkAm^>  = 

In  detail,  the  mean  of  the  estimated  centroid  is 


<<=>  = a( 


1 + H >wn2  ) -Z  £k 

k,-t  k l 


(id 


and  its  covariance  matrix  is 


<AS*S>  = / ] fcvEr.  - (sr,  + rv  J + ss]^,  /N2  . 

k,£ 


(12) 


For  the  second  moment,  the  mean  and  the  fourth-rank  covariance 
tensor  are 


<i>  + ) ~ 'y  ]~k~k  ^ Vlct ) > (13) 


k ,1 


<^l>  [£y£v£f.£f.  - + £>£>”)  /n2  , (14) 

k ,1 


in  which  juxtaposed  vectors  and  matrices  denote  outer  products 
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The  diagonal  elements  of  the  detection  matrix  p,  are  the 
variances  of  the  number  of  photocounts  at  the  detector  out- 
puts. Some  properties  of  the  variances  were  discussed  in 
Ref.  5,  where  Mandel's  semiclassical  theory  of  photodetection1 3 
was  applied  for  several  types  of  laser  illumination  and  target 
surface  roughness.  The  off-diagonal  elements  of  p,  involve 
correlations  of  photoelectrons  among  pairs  of  detectors  in  the 
array.  This  correlation  was  demonstrated  first  by  Hanbury 
Brown  and  Twiss.14 

The  elements  of  p.  involve  fourth-order  correlations  of 
the  optical  field,  corresponding  to  second-order  correlations 
of  the  irradiance,  and  are  therefore  expressible  as15 


^ t 6k^ + <nk^<n-*,^>  su 


(15) 


where  6^  is  the  Kronecker  delta  and  is  an  integrated 

cross  correlation  of  the  instantaneous  intensity  l(r,t)  at 
the  k-th  and  -t-th  detectors,  at  different  times  t,t' : 


<&.  =<AIkAIt>/^<TT>  • <16> 

The  overbars  denote  spatio-temporal  averaging  over  the  detector 
apertures  and  integration  time. 

The  first  term  in  (15)  represents  the  quantum  noise  in  the 
count  of  photons  from  the  image  and  background  radiation  and 
the  detector  dark  current.  The  other  term  reflects  random  fluc- 
tuations in  the  energy  W and  involves  the  covariance  of  the 
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energy  at  pairs  of  detectors.  In  Ref.  5,  was  defined  as 

the  contrast  of  the  measured  fading  of  the  fluctuating  irra- 
diance,  generalizing  the  quantity  known  as  contrast  of  the 
laser-speckle  pattern. 

For  optical  fields  that  can  be  modeled  as  complex  Gaussian, 
weakly  stationary  processes,  is  related  to  their  spatio- 

temporal  degree  of  coherence  by 

CjU  - i-Wt-OI’  • (17) 

For  monochromatic  illumination  and  rough  targets,  for  which 
fluctuations  in  W are  strongest,  the  image  field  is  such  a 
process,  over  space.  The  degree  of  coherence  is  then  only  spa- 
tial and  is  fully  determined  by  the  point  spread  function  of 
the  telescope.5  For  the  same  illumination  but  targets  with 
smooth  surfaces,  there  are  no  random  fluctuations  in  energy 
and  the  last  terrain  (15)  is  absent.  This  term  is  also  zero 
for  temporally  noncoherent  illumination  with  coherence  time 
much  shorter  than  the  integration  time. 
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DISCUSSION 


To  aid  in  interpreting  the  general  expressions  for  the 
performance  and  accuracy  of  the  estimators,  it  is  useful  to 
relate  them  to  the  SNR  of  the  array  as  a whole.  This  compo- 
site sensor  registers  a random  total  reduced  count  ^Tn^, 
with  mean  N and  variance  XI  <I^le  si9nal“to-noise  ratio 

at  the  output  of  the  array  is 


R = 


(18 


The  fractional  covariances,  summing  to  unity,  are 


Tc  l 


= “w/2  “ij 


(19 


The  correlation  between  an  individual  reduced  count  n^  and 
the  total  count,  normalized  to  the  overall  variance,  is 

vk  = ^"vAk£  • (2l 

l 

The  centroid  of  this  fluctuation  count  is  located  at 

3 (2: 

k k,£ 

and  its  second  moment  is 
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In  terms  of  these  statistically  meaningful  quantities. 


the  mean  and  accuracy  of  the  estimators  are  expressible  by 


<£>  = l + (5  “S)/r2  > 


(23) 


<*8&s>  -[2  £kV- 


~ qq  + (s 


-q)(s  -q)l 

w <v  ~ I 


/R2  ; (24) 


<i>  = M + (M-Q)/R2  , 


(25) 


<A^i>  =C^  Ek^kAk^£t-SQ+  (M-Q)(M-Q)1/R2  . (26; 

^ -U  P J 


The  rms  errors  in  the  estimates  are  inversely  proportional  to 


the  overall  SNR;  the  biases  are  inversely  proportional  to  the 


square  of  the  SNR. 


Equations  (23)  and  (24)  can  be  applied  to  a tracking  sys- 
tem in  which  the  displacement  of  the  image  centroid  from  the 
optical  axis  is  used  to  control  the  monopulse  laser  radar  track- 
ing the  TAP.  The  system  can  be  conceptually  similar  to  typical 
ones  employed  for  tracking  point  targets,  in  which  the  peak  po- 
sition of  the  imaging  system's  point  spread  function  is  tracked 
with  a four-quadrant  detector.4  in  what  follows,  these  equa- 
tions are  specialized  to  cases  in  which  the  detection  matrix  (x 


retains  only  its  diagonal  terms. 


■ . L 


RIVERSIDE  RESEARCH  INSTITUTE 


where  sw  is  the  centroid  of  the  energy  fluctuations 
<nk>S|>i*|S  • In  the  absence  of  the  photon-bunching  effect, 
only  the  first  term  in  (29)  remains,  and  this  bias  becomes  neg- 


ligible when  the  signal  count  N becomes  substantial.  Even 


when  strong  energy  fluctuations  are  present,  it  should  be  noted 


that  s - sw  is  a measure  of  the  nonuniformity  of  the  image 
irradiance  and  K2  is  generally  less  than  unity.  From  (28), 


one  observes  that  K-1  can  be  interpreted  as  the  SNR  at  the 


array  output,  when  it  is  limited  by  the  laser-speckle  noise. 


This  SNR  increases  as  the  square  root  of  the  number  of  detectors 


sensing  the  image  intensity;  for  Nyquist  sampling,  this  equals 


the  number  of  degrees  of  freedom  in  the  image. 


The  diagonal  detection  matrix  also  yields  for  the  covar- 


iance matrix  of  the  centroid  estimate 


<AcAc>  = (M— J5s)/N+  (Ma  + ss)(Nn/N2)+K2  [(^-8^)+  (s  - s^  (s  - £w)] 


(30) 


in  which  M - ss  is  the  image's  intrinsic  moment  of  inertia, 
referred  to  its  own  centroid,  ^ is  the  corresponding 

moment  of  the  speckle  noise,  and  1^  that  for  the  array.  The 
geometrical  factors  depend  on  image  size  orientation,  and  loca- 


tion; they  do  not  change  with  laser  power,  provided  the  illumi- 


nation is  uniform.  The  coefficients  are  the  inverse  — square 


SNR  of  the  sensor,  when  limited  by  the  signal  3hot  noise,  back- 


ground and  dark  current  noise,  and  laser  speckle.  The  last  of 


j./, 
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these  is  insensitive  to  changes  in  the  laser  power,  while  the 
others  can  be  reduced  by  increasing  the  power.  For  large  sig- 
nals, the  laser-speckle  noise  limits  the  system’s  accuracy. 

For  fixed  longitudinal  geometry  of  the  optical  system  and  given 
laser  coherence,  this  limitation  can  be  eased  by  increasing 
the  space -bandwidth  product  of  the  imaging  system  by  an  increase 
of  the  telescope  aperture.  This  not  only  improves  the  diffrac- 
tion-limited resolution  of  the  imaging  system,  but  also  produces 
more  degrees  of  freedom  in  the  image,  thereby  increasing  the 
speckle-noise-limited  SNR.5  For  nonfluctuating  targets,  the 
last  term  in  (31)  is  absent. 

The  square  root  of  the  trace  of  the  matrix  <(AcAc]>  , is 
the  rms  error  in  the  centroid  estimate.  Assuming,  for  simplic- 
ity, a uniformly  illuminated  and  reflecting  target  and  cross- 
spectrally  pure  light,  one  has  aw  = s and  MK,=M.  The  accu- 

'"wW  ^ ^ 

racy  of  the  centroid  estimate  is  then  given  by 

o/rz  = [K2  +1/N+  (NN/fcT2)(rj[  + a2)/r\y/2  , (31) 

whers  a is  the  standard  deviation  of  the  centroid  location, 
rj  and  rA  are  the  optical  radii  of  gyration  of  the  image  and 
array,  respectively,  and  s is  the  distance  of  the  image  cen- 
troid from  the  optical  axis.  Figure  2 is  an  asymptotic  log-log 
plot  of  normalized  rms  error  as  a function  of  signal  N.  It 
shows  the  combined  effects  of  total  signal  photocounts,  back- 
ground and  dark  current  counts,  and  laser-speckle  noise,  as 


*V  V -V  . 
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well  aa  the  impact  of  geometric  factors,  including  image  size, 
centroid  displacement,  and  sensor  size.  As  the  signal  increases 
from  weak  levels,  the  accuracy  is  limited  first  by  background 
light  and  dark  current,  and  varies  with  image  location.  At 
intermediate  signal  levels,  the  image  size  and  the  shot  noise 
determine  the  accuracy.  For  strong  signals,  the  accuracy  be- 
comes independent  of  signal  level,  being  limited  by  speckle 
noise,  but  still  depends  on  the  size  of  the  image. 


_V  I 
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CONCLUSIONS 

The  angular  position,  shape,  and  orientation  of  extended 
targets  can  be  estimated  by  use  of  monopulse  laser  radars,  with 
target  images  sensed  by  a photon  counting  array.  The  processing 
involves  conversion  of  the  output  counts  into  measurements  of 
the  image  centroid  and  optical  moment  of  inertia. 

A Taylor  series  expansion,  valid  for  small  fluctuations 
of  the  photocounts  about  their  means,  was  used  to  derive  the 
means  and  covariance  matrices  of  the  estimators, to  the  second 
order  of  smallness.  It  has  been  found  that: 

a.  Both  estimators  are  biased;  the  biases  are  inversely 
proportional  to  the  square  of  the  overall  SNR  at  the  sensor 
output. 

b.  The  rms  errors  are  inversely  proportional  to  the  SNR. 

c.  The  biases  and  rms  errors  of  the  estimators  depend 
on  the  detection  matrix,  comprised  of  the  covariances  between 
the  photocounts  at  the  outputs  of  pairs  of  detectors  in  the 
array.  The  detection  matrix  becomes  diagonal,  giving  the  var- 
iances of  the  photocounts  at  the  individual  detector  outputs, 
for  either  non fluctuating  targets  or  noncoherent  illumination, 
or  also  for  fluctuating  targets  generating  Gaussian  image  fields, 
if  the  detectors  sample  them  at  the  Nyquist  rate. 
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d.  For  diagonal  detection  matrices  and  nonfluctuating 
targets , the  bias  and  rms  error  of  the  centroid  estimator  de- 
crease with  increasing  laser  power;  for  fluctuating  targets 
they  cannot  be  decreased  below  the  limit  imposed  by  laser- 
speckle  noise.  This  limitation  improves  with  the  number  of 


degrees  of  freedom  in  the  image , which  depends  on  the  diffrac 


tion-limited  resolution  of  the  imaging  system. 
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Fig.  1. 


Figure  Captions 


Image  Plane.  I = image,  = position  of  k-th  detec- 
tor, = number  of  photocounts  at  r^,  c = cetroid  esti- 
mator,  <c;>  = mean  value  of  centroid  estimator, 

Ac  = c - <(c^>,  and  s = true  image  centroid. 


Fig.  2.  Normalized  rms  error  of  c vs  mean  signal  counts; 


asymptotic  approximations.  The  steepest  asymptote 
is  background-noise -limited:  px  = r^fNjjfr^  + s2 

Nx  = ( rj^  + s 2 )Njj/r ^ . The  middle  asymptote  is  for  sig- 


nal-shot-noise-limited  operation:  p = K,  N = K 2. 


The  horizontal  asymptote  represents  laser-speckle- 


li.xLted  operation. 


! ■*  *.  ».  » * . * • • * ■ " » 

. • . *. * . ^ _ « • . • . • _ 

V 


- ” » • • - __  _ 

- / v v v • 


If 


SECURITY  CL  ASSIF  1C  AT  ION  Of  TmS  PACC  flO«.(i  flolt  r«>n« 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


MEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  OOVT  ACCESSION  NO 


2.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (end  Submit) 

Photocounting  Image  Tracking  of 
Fluctuating  Targets 


S.  TYPE  OP  REPORT  • PCRIOO  COVERED 

Technical  Report 


f.  PERPORMINO  ORO.  REPORT  NUMRI 

Tl/364-3-65 

S.  CONTRACt  OR  GRANT  NUMBER!*) 


T.  AUTHOR!*) 

Marek  Elbaum 


DAAK40-76-C-0500 


•0.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  * WORK  UNIT  NUMRERS 

ARPA  Order  2281,  Amd.l^ 
Tech.  RQMT  6004 
Par.  3.4 


S.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Riverside  Research  Institute 

80  West  End  Avenue 

New  York,  N.Y.  10023 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

MIRADCOM 

Redstone  Arsenal,  Alabama  35809 


12.  REPORT  OATE 

February 


1977 


IS.  NUMBER  OF  PAGES 

116  + vi 


14.  MONITORING  AGENCY  NAME  a AOORESS (II  dill* ml  from  Ctnlltlllni  Olllet) 


■ S.  SECURITY  CLASS,  (of  Ml*  r.port) 


UNCLASSIFIED 


ISa.  OECLASSI  FI  CATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  I hit  RtptU) 


IT.  DISTRIBUTION  STATEMENT  <tl  Iht  abairtcl  tnltrtd  In  Block  JO.  II  dllltitnl  from  Rip on) 


Unlimited  Distribution 


v«.  supplement  ary  notes 


19.  KEY  WORDS  ( Conllnu • on  tavraa  • aida  (/  n«(M<«ry  and  IdantUy  by  block  numbor) 

Laser  Radar  Laser  Speckle 

Position  Estimation  Photocounting 

Tracking  Noncoherent  Detection 


20.  ABSTRACT  (Continue  on  reetrte  tidt  II  nectteety  tnd  Identity  by  block  numbet) 

A theory  of  estimation  of  angular  position  and  other  attribute 
of  optically  rough  and  smooth  targets  with  a monopulse  laser 
radar  is  developed.  It  applies  to  systems  deriving  informa- 
tion about  the  target  position  by  sensing  its  image  with  an 
array  of  noncoherent  detectors.  The  theory  develops  quanti- 
tative formulations  of  the  fundamental  limitations  imposed 
upon  measurement  accuracy  by  the  shot  noise  arising  from  both 


wcu»mr  cimvhohok  or  ihh  Mtinwnn  i>*»  «■!««»<) 


StCUMTV  CLASSIFICATION  OF  THIS  Of  !"•«« 


y <• 


